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ABSTRACT 


Passive bistatic radar (PBR) and passive multistatic radar (PMR) use opportunistic 
transmitters to detect and locate targets. In this thesis, a maritime scenario was modeled 
with merchant vessels serving as multiple non-cooperative opportunistic transmitters, 
while a frigate warship equipped with Electronic Warfare (EW) and Direction 
Finding (DF) receivers takes on the role of the receiver in a PBR/PMR configuration. The 
targets are assumed to be the generic Formidable-class frigate. 

A MATLAB model is developed to simulate the operating environment and 
passive detection and location process. Detection coverage is investigated to propose 
optimal PBR/PMR configurations and geometry, while elliptical and hyperbolic target 
location methods are explored to quantify the effects of PBR/PMR parameters and 
geometry on target estimated location uncertainty. 
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I. INTRODUCTION 


A. OVERVIEW 

Since the concept of radar engineering was first demonstrated in 1904, progress in 
radar technology has been driven by growing requirements for radar performance and 
rapidly changing operating environment. A long detection range and wide coverage, 
measurement accuracy, greater system capacity, and an ability to operate with the 
presence of interference are some of the fundamental radar requirements that have been 
established over the past few decades [1]. These radar performance characteristics— 
together with the need to detect, separate, classify, locate, and track sources of emissions 
in multi-target environments—triggered the development of passive radar detection and 
location techniques. 

The vast majority of today’s deployed radar systems are monostatic, that is, the 
transmitting and receiving antennas are collocated. Despite the advancements in 
antennas, transmitters, receivers, and processing technology, as well as passive radar 
systems, conventional monostatic radar remains a double-edged sword—whereby it 
detects targets, but radar transmission makes it vulnerable to detection—and may not be 
the best option to address certain operational scenarios. A promising solution is to use 
multiple radar transmitting and receiving sites to exploit spatial advantage for 
coordinated target detection. Multi-site radars can be broadly classified into bistatic radar 
and multistatic radar. Bistatic radar is a radar system where the transmitter and receiver 
are located at different sites [2], Similarly, a multistatic radar system utilizes multiple 
spatially separated transmitter and receiver sites where the target information from all 
receivers is fused [1]. Passive bistatic or multistatic radar capitalizes on transmitters of 
opportunity to detect and locate sources of transmission or targets without deliberate 
emissions. The illuminators are not limited to radar signals and include (but are not 
limited to) analog TV, FM radio, digital video broadcasting terrestrial (DVB-T), digital 
audio broadcasting (DAB), cellular network, WiFi, and Global Positioning System (GPS) 
satellite signals [3]. 
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In naval operations, targets employing low radar cross section (RCS) and radar- 
absorbent material (RAM) design methods, coupled with the use of highly sensitive 
electronic warfare (EW) receivers on warships, have changed the nature of the game. 
Military ships are pressured to limit transmissions to avoid detection by highly sensitive 
EW receivers. Furthermore, the use of low probability of intercept (LPI) radar in a 
monostatic configuration results in weak returns from low-RCS targets, restricting 
detection capability and compromising situational awareness. This thesis explores the use 
of multiple pairs of passive bistatic radar (PBR) to detect low-RCS targets using 
opportunistic transmissions as a possible solution and to validate the findings in [4]. 
Figure 1 shows the bistatic radar geometry for a PBR pair. The direct line-of-sight (LOS) 
distance between the transmitter (Tx) and receiver (Rx) is known as the baseline and is 
denoted by L . The distance between the transmitter and target is denoted as R T while 

the distance between the receiver and the target is denoted as R r . The bistatic angle /3 is 
the angle subtended between transmitter, target, and receiver. 

The proposed PBR approach offers potential advantage in the detection of 
stealthy, low-RCS targets which are designed to minimize monostatic radar echoes. RCS 
returns of stealthy ships vary with bistatic angle />' and can be sufficiently large at certain 
return directions. PBR capitalizes on this characteristic using opportunistic transmissions 
to detect low-RCS targets. Being a passive system, PBR allows the receiver to remain 
covert, making it more resilient to detection and electronic attack (EA)—in the form of 
jamming and anti-radiation missiles (ARMs). The ability to leverage available 
transmission and to detect targets passively also serves to enhance situational awareness, 
thereby advancing one’s position in the observe-orient-decide-act (OODA) loop during 
an operation [4], It is also advantageous to use multiple transmitters at separate locations 
for detection as it adds spatial diversity, which enhances detection accuracy and aids in 
removing clutter, interference, and potential system errors. Lastly, the PBR system 
proposed requires no additional equipment as all necessary hardware is currently 
available on most naval ships. 
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Target (Tgt) 



Figure 1. Bistatic radar geometry. 


B. HISTORY 

The concept of Radio Detection and Ranging (RADAR) was first demonstrated in 
1904 by German inventor Christian Hiilsmeyer in a monostatic setup [5]. Following this, 
radar experiments in the United States, the United Kingdom, France, Italy, Russia, and 
Japan were carried out predominantly with bistatic radar operating in the forward¬ 
scattering configuration [6]. However, the invention of the radar duplexer in 1936 
addressed the issue of transmitter-receiver isolation and broadened the application of 
single-site radar on aircraft, ships, and mobile ground units. By the end of World War II, 
bistatic radar was displaced by monostatic radar, with the former experiencing periodic 
resurgences [5]. 

1. First Resurgence 

The 1950s saw renewed interest in bistatic radar with developments in missile and 

satellite detection, location and tracking, semi-active homing missiles, hitchhiking, and 

second-generation forward-scattering fences [2], During this period, advancement in 
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radar theory led to a better appreciation of concepts on match filtering, ambiguity 
functions, statistical theories on detection, moving target indicator (MTI) radar, and 
synthetic aperture radar (SAR) [7-9]. The term bistatic radar originated with K. M. 
Siegel and R. E. Machol in 1952 [10]. 

2. Second Resurgence 

The development of counter-measures to anti-radiation missiles (ARMs) and 
emitter location-directed artillery in the 1970s resulted in the second resurgence in 
bistatic radar. With a dual or multiple site configuration, the effectiveness of electronic 
support measures (ESM) directed attacks can be reduced by locating the transmitter away 
from the receiver or into a less vulnerable sanctuary [5]. The advent of digital processing 
techniques during this period increased the processing capability in MTI operation modes 
and allowed real-time airborne SAR mapping [7], 

3. Third Resurgence 

Research on bistatic space-time adaptive processing (STAP) to address moving 
clutter and concepts to improve bistatic SAR images signaled the start of the third 
resurgence. It was also during this period that passive bistatic radar surfaced as a possible 
counter stealth technique. The idea for PBR is to utilize commercial broadcast signals in 
bistatic or multistatic configurations to detect low-RCS targets [5]. Since then, several 
PBR systems have been developed and evaluated for air surveillance. Some notable PBR 
air surveillance systems are listed in Table 1. 
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Table 1. Parameters of significant passive bistatic radar programs designed 
and tested for air surveillance.* After [5]. 


System 

Silent 

Sentry™ 

TV-Based 
Bistatic 
Radar (I) 

TV-Based 
Bistatic 
Radar (II) 

FM Radio - 
Based 
Bistatic 
Radar 

Multistatic 

HDTV-Based 

Radar 

Developer 

IBM, now 
Lockheed 
Martin 

Univ College 
London 

DERA, United 
Kingdom 

NATO 

SAIC, U.S. 
Army 

Decade 

configuration 

1980-2000 

1980 

1990 

2000 

2000 

Transmitter 

operation 

Multistatic: 

Rx: 1 

Tx: up to 6 

Bistatic 

Bistatic 

Near-forward 

scatter 

Bistatic 

Multistatic: 

Rx: 4 

Tx: 1 

Baseline 

100 km 
typical 

12 km 

150 km 

50 km 

10 km 
typical 

Target 

Aircraft 

Missile 

launches 

Aircraft 

Aircraft 

Aircraft 

Aircraft 
below ~ 5000 
ft 

Target data 

Range 

Doppler 

Bearing 

Range 

Bearing 

Doppler 

Bearing 

Range 

Doppler 

Bearing 

Range 

Doppler 

Measured 

performance 

R m = 100 km 

150 km 

2- D tracks on 

A/C 

3- D tracks on 

missile 

launches 

R m ~ 25 km 
Occasional 
A/C 

detections, but 
mostly 
negative 

R m ~ 160 km 
Detections on 
high and 
medium 
altitude A/C 
but only 1/3 
tracked 

R m ~ 175 km 
Achieved with 
innovative 
direct path 
excision 

R m ~ 30 km 
Target location 
via multi- 
lateration 
Ghost excision 
via Doppler 
association 

Status 

Version 3 for 
sale to U.S. 
Government 
for < 1 million 
dollars 

Work 

continues 

Program 

ended 

Program 

ended 

Program 
continuing 
possible in a 
multistatic 
mode 

Test phase 
complete 
Awaiting 
evaluation/ 
funding 


is the equivalent maximum monostatic range defined as (WZ =K , where R r is the 


receiver-to-target range and R T is the transmitter-to-target range. 


5 







C. RECENT DEVELOPMENTS IN BISTATIC AND MULTISTATIC RADAR 

Despite research progress on improving the detection, classification, and location 
performance of passive bistatic radar systems, their reliance on transmissions of 
opportunity and the restricted geometry has limited their application. This fundamental 
requirement continues to stimulate research and experimentation on opportunistic 
illuminating sources, their optimum configuration, applicability, and performance in 
various operational scenarios. As with all radar systems, improvement in interference and 
clutter rejection, target detection, classification, and location and tracking accuracy are 
current areas of interest in the field of passive bistatic radar and passive multistatic radar 
(PMR). The recent research topics in bistatic and multistatic radar can be classified as 
bistatic and multistatic system configuration, forward-scattering radar, and multiple-input 
multiple-output (MIMO) radar. 

When the concept of bistatic and multistatic radar was first introduced, dedicated 
radar transmitters were used as transmission sources [11] before transmitters of 
opportunity were employed as illuminating sources. Transmitters of opportunity in the 
very high frequency (VHF) and ultra-high frequency (UHF) band, such as FM broadcast, 
TV broadcast, DAB, DVB-T, and cellular network signals continue to be common 
illuminating sources used in detecting airborne, land, and maritime targets [11-19]. 
Conversely, studies on the use of high frequency (HF) band signals as opportunistic 
transmitters only started recently. HF signals present advantages of long range detection 
and coverage, propagation beyond the radar horizon, and improved detection of stealth 
targets, which enhances the PBR/PMR’s early warning capability [20, 21]. A list of 
common transmission sources and their typical parameters are given in Table 2. 

Apart from terrestrial sources of transmission, there has been recent interest in 
using satellite transmissions to detect airborne targets. The SABER-DEMO platform 
recently demonstrated its ability to detect aircraft passively using signal processing 
techniques to process weak satellite sources of transmission [22], 
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Table 2. Signal parameters for typical passive radar illumination sources. 

From [23], 


Transmission 

Frequency 

Modulation, 

Bandwidth 

P,G t 

Power Density (Wm‘) 

®= p,Gt 

4^Tj 2 

HF broadcast 

10-30 MHz* 

DSB AM, 9 kHz 

50 MW 

-67 to -53 dBW m 2 atri 
= 1000 km 

VHFFM 

(analogue) 

-100 MHz 

FM, 50 kHz 

250 kW 

-57 dBW nf 2 at ri = 100 
km 

UHFTV 

(analogue) 

-550 MHz 

Vestigial-Sideband AM 
(vision); FM(sound), 

5.5 MHz 

1 MW 

-51 dBWm- 2 atr 1= 100 
km 

Digital audio 
broadcast 

-220 MHz 

digital, OFDM 220 
kHz 

10 kW 

-71 dBWm- 2 atr 1= 100 
km 

Digital TV 

-750 MHz 

digital, 6 MHz 

8 kW 

-72 dBW m' 2 at n = 100 
km 

Cellphone base 
station (GSM) 

900 MHz, 

1.8 GHz 

GMSK, 

FDM/TDM A/FDD 200 
kHz 

100 W 

-81 dBW m' 2 at r! = 10 
0km 

Cellphone base 
station (3G) 

2 GHz 

CDMA 5 MHz 

100 W 

-81 dBW m 2 at r! = 100 
km 


♦Appropriate frequency will depend on time of day. 


Another topic of interest is the use of inverse synthetic aperture radar (ISAR) 
processing techniques with forward-scattering radar (FSR) for target detection and 
parameter extraction. FSR is the earliest form of bistatic radar where target detection 
occurs at the transmitter to receiver baseline. The target’s radar cross section is enhanced 
in the forward-scattering configuration due to Babinet’s principle [2], Using the constant 
false alarm rate (CFAR) approach, target detection and parameter extraction can be 
accomplished in real time [24], An experiment conducted at Forsvarets forskningsintitutt 
(FFI) explores the difference in ISAR ship signatures in the forward and back-scattering 
configurations. The results from the study show that forward-scattering returns produce 
more accurate ISAR signatures than back-scattering returns as the RCS for forward¬ 
scattering is usually stronger than back-scattering. The difference between forward and 
back-scattering ISAR signatures can be fused to improve ship identification and 
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classification [25], Apart from studies on FSR processing techniques to improve target 
detection and extraction, considerable research has been conducted on ultra-wideband 
(UWB) FSR for its ability to reduce sea clutter in maritime applications [26, 27]. 

Recent advances and interest in MIMO radar systems can be attributed to its 
potential for detection and location of targets in bistatic or multistatic configurations. 
Spatial diversity in MIMO radar systems exploits the differences in target cross section in 
detecting and extracting target parameters such as angle of arrival (AOA) and Doppler 
frequencies [28]. Furthermore, coherent processing in MIMO systems improves target 
location accuracy [28, 29]. 

As with all studies in the radar domain, current bistatic and multistatic research 
areas are motivated by the need to improve detection capability and measurement 
accuracy while reducing or mitigating the effects of unwanted interference. 

D. THESIS OBJECTIVE 

The primary objectives of this thesis are as follows: 

1. Generate a MATLAB model that computes a system’s detection 
performance given the target’s path. 

2. Examine low-RCS maritime target detection coverage and performance 
using multiple pairs of bi static radar. 

3. Investigate low-RCS maritime target location accuracy using elliptical and 
hyperbolic target location methods. 

The EM simulation software FEKO will be employed to model and analyze the 
RCS scattering properties of the maritime target while MATLAB will be used to model 
and simulate the operating environment and PBR target detection and location. The 
MATLAB model is validated against open source literature by using simple targets with 
known RCS return characteristics prior to using the models to generate simulation results. 

The MATLAB model generated for objective (1) requires the user to provide an 
Excelfile with the target’s position and velocity vector components at each time step, 
after which detection performance parameters will be computed. This model will be a 
fully automated MATLAB model that calls FEKO to compute precise RCS returns given 
the exact incident and receive angles. Detection coverage results from objective (2) will 
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be compared against findings in [4] and used to propose the optimal PBR configuration 
and geometry for maximum detection coverage. Target location estimation results from 
objective (3) will be used to examine the effects of PBR geometry on location error 
distribution. 

The maritime scenario with a single receiver and multiple transmitters depicted in 
Figure 2 applies to all models. The problem setup assumes a warship (receiver) equipped 
with broadband EW receiver and direction finding (DF) capability deployed to monitor 
maritime traffic flow in the Straits of Singapore. The targets are assumed to be low-RCS 
targets with infrequent transmissions and/or operating with LPI radars, while civil marine 
radars (CMR) on merchant ships will serve as opportunistic illuminators in a PBR 
configuration [4], Regulation 19 of SOLAS Chapter V requires all merchant and 
warships to carry automatic identification systems (AISs), which share information on the 
ship’s identity, position, course, speed, navigation status, and safety-related information 
[30], 



-> 

-> 


Merchant Ships (Tx) 

Maritime AIS and Nav Radar Tx (direct) 
Nav Radar Tx (indirect) 


Figure 2. Pictorial representation of the PBR setup. After [4], 
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E. THESIS CHAPTER OUTLINE 


Chapter I introduced the concept, history, and recent developments in passive 
bistatic radar. The goals and end products of the thesis are also detailed here. 

Chapter II provides the theoretical background on PBR. The corresponding 
parameters employed to develop the necessary MATLAB model are introduced. The 
bistatic range equations, detection contours, and mathematical concepts related to target 
location and uncertainty are covered. 

Chapter III presents the design approach in modeling the scenario and the 
problem setup using FEKO and MATLAB. The results from verifying the detection 
coverage models using simple targets and PBR geometry against known results are 
covered in detail. The methodology used to verify target location model is also discussed. 

Chapter IV uses the model generated to examine detection coverage for different 
PBR transmitter-target-receiver geometries. Simulation results are compared against 
findings in [4] and used to propose PBR configurations and geometry for optimal 
detection coverage. The results for elliptical and hyperbolic target location methods and 
their corresponding uncertainty ellipse are presented. The effects of transmitter-target- 
receiver geometry target location accuracy are also explored. 

Chapter V summarizes the research findings and suggests further work to improve 
the models and multiple PBR detection, parametric extraction, and location capability 
and accuracy. 
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II. BISTATIC RADAR THEORY 


A. DEFINITION 

Bistatic radar refers to a radar system where the transmitter and receiver are at 
sufficiently different locations such that the angles or ranges from those locations to the 
target are significantly different [31]. The basic bistatic configuration and parameters are 
defined in Figure 1. 


B. RANGE EQUATION 

The bistatic radar range equation (BRRE) gives the received power at Rx as a 
function of the system parameter, target scattering properties, and engagement geometry. 
Solving the BRRE for the range product gives [2] 


where 


k = (R T R R ) max 


P t G t G r A 2 g B F R F R G p 
(AxfkT'B'iS/W^LrLt 


( 1 ) 


R t = transmitter-to-target range, 

R r = receiver-to-target range, 

P T = transmitter power output, 

G t = transmitting antenna power gain, 

G r = receiving antenna power gain, 

A = wavelength, 
c j b = bistatic target cross section, 

F t = pattern propagation factor for transmitter-to-target path, 
F r = pattern propagation factor for target-to-receiver path, 

G p = processing gain, 
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k = Boltzmann’s constant (1.38x10 23 J/K), 

T s = receiving system noise temperature, 

B n = noise bandwidth of receiver’s pre-detection filter, sufficient to pass all 
spectral components of the transmitted signal, 

(,S / /V) mm = signal-to-noise power ratio required for detection, 

L r = transmitting system losses (>1) not included in other parameters, 

L r = receiving system losses (>1) not included in other parameters, 

K = bistatic maximum range product. 

In the bistatic range equation, the maximum range product RpR R replaces R 2 in 
the monostatic range equation where R = Rj = R r is the monostatic transmitter-to-target 

and target-to-receiver range. The difference between the transmission path and receiving 
path results in significant differences between monostatic and bistatic radar operation. 

One of the differences is that monostatic contours of equal signal strength are 
constant range circles, while detection contours for bistatic radar are defined by ovals of 
Cassini. An oval of Cassini is defined as a locus of points where the product of the 
distance from two fixed points is constant. Figure 3 shows the Cassini oval for two fixed 
points ( F 1 and F 2 ) separated by a distance of 2 c. 


Applying the concept of Cassini ovals to the bistatic triangle in Figure 4 with 
baseline L and range product R t R r , an expression for constant signal-to-noise (SNR) 
power ratio can be derived by writing Eq. 1 as [2] 


S/N = 


K 

p 2 r)2 

k t k r 


( 2 ) 


where SIN is the signal-to-noise power ratio at R r and R r , and K is the bistatic radar 
constant 
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From the geometry in Figure 4, R T and R r are converted to polar coordinates 


(, r,9 ) using the law of cosines: 


R; = 


(r 2 +L/^) + rLcos0, 


R 2 r = (^ 2 + ~ rL cos 6, 


(4) 

(5) 


where the origin is at the midpoint of the baseline. Substituting Eq. 3, Eq.4, and Eq. 5 
into Eq. 2 gives an expression that defines constant SNR contours [2]: 


S / N = 


_ K _ 

(r 2 +L 2 /4) 2 - r 2 l} cos 2 6 


( 6 ) 


Signal-to-noise ratio contours generated using Eq. 6 for 10 dB <S/N< 30 dB and 
K = 30L 4 are given in Figure 5. 



Figure 5. Ovals of Cassini, contours of constant SNR (dB), with K = 30L 4 . 

After [2], 
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Given that L , ( R t +R r ) and 0 R are obtained and measured from the 


opportunistic transmitter and receivers, R T and R r are computed as [2] 


R t = (R 1 - +L 2 + 2 R r L sin 0 R ) 1/2 , 

R (/? r +^) 2 ~L 2 
2(Rj + R r + L sin 0 R ) 


(7) 

( 8 ) 


Using the law of cosines on the bistatic triangle in Figure 4 yields 


P = cos 


-i 


R 2 +R 2 r 


r2\ 


2 R t R r 


(9) 


V J 

When L = 2 yfic , the oval forms a lemniscate with cusp at the origin. The ovals of 
Cassini in Figure 5 define three operating regions for bistatic radar: 


1 . L > 2\fic with Rj. » R r . Receiver centered region. 

2. L > 2\[k with R r « R,. Transmitter centered region. 

3. L < 2\[k . Cosite region or receiver-transmitter-centered region. 

In cases where the target echo signal strength is weak, non-coherent pulse 
integration performed after the envelop detector increases SNR by a factor of ^//V where 
N is the number of pulses integrated. Improvement in SNR by pulse integration is a 
form of processing gain. The number of pulses integrated over a period of t seconds is 
calculated as 


where 


N = PRFxTOTxt 
60 


TOT = 


360 co_. 


( 10 ) 

( 11 ) 


6 mb = 3 dB azimuth beamwidth, 


co scan = scan rate (rpm), 

TOT = time on target, 

PRF = pulse repetition frequency, 

t = integration period. 


15 



For the purpose of this thesis, the minimum difference between the noise floor and signal 
level is assumed to be 10 dB, that is, minimum SNR is 10 dB. 

C. TARGET LOCATION EQUATIONS 
1. Bistatic Radar Trigonometry 

The AOA of the target echo signal 0 R and target-to-receiver range R r are 

required to define the target’s location with respect to the receiver in a PBR 
configuration. The AOA of the echo signal can be measured directly; however, the target- 
to-receiver range cannot be measured directly and needs to be calculated by solving the 
parameters of the bistatic triangle (Figure 4). 

To solve for target-to-receiver range R r and the rest of the bistatic triangle 
parameters requires measuring and knowledge of the following: 

• Baseline range from transmitter position(s) and receiver position, L , 

• AOA of target echo signal at the receiver, 0 R , 

• Transmitter-to-target and target-to-receiver range sum, (R t +R r ). 

The range sum ( R T +R r ) can be estimated using the direct and indirect method as 
illustrated in Figure 6. In the direct method, the receiver measures the time delay A T rt 
between the reception of the transmitted pulse and the target echo. The range sum can 
then be expressed as a function of the time delay A T rt and the baseline range [2]: 

(R t + R r ) = cAT rl + L. ( 12 ) 

In the indirect method, the receiver measures the time delay AT tt between the 
transmission of the pulse and the reception of the target echo. The range sum in this case 
is a function of the time delay AT tt [2]: 

(Rj. +R r ) = cAT tt . (12) 

The indirect method requires receiver and transmitter clocks to be synchronized while the 
direct method can be used with any transmitter configurations given LOS between the 
trans mi tter and receiver. 
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Referring to Figure 7, the direct and reflected pulses received must be resolvable 
such that 


t 


reflected 


>f direct +T 


(14) 


a) Timing Sequence 


TIME 


Pulse emitted 
by Tx 


Pulse arrives Pulse arrives 

at Rx at Tgt 


Tgt echo 
arrives at Rx 



b) Direct Method 


Ry -I - ~ f 


A 1 . 


Start Clock 


Stop Clock 


c) Indirect Method 


R 


+ Rk-cA7 


Start Clock 


AT, 


Stop Clock 


Figure 6. Timing sequence diagram for direct and indirect method for 
calculating range sum (R T +R r ) . From [2], 

17 
























Figure 7. Reception of direct and reflected pulses. 


To derive the mathematical relationship between the bistatic triangle parameters, 
first consider the elliptical iso-range contours on a bistatic plane such that each concentric 
ellipse is determined by 

Ri R r 2a, (15) 

where a is the semi-major axis length of the ellipse (Figure 3). The eccentricity of the 
ellipse e is therefore defined as 


e = ■ 


2 a 
L 


e = 


(■ Rt + Rr) 


(16) 


(17) 


Given the measurement of ( R T + R h ,), L, 0 R and using the law of cosines on the 


bistatic triangle in Figure 4, 

= R 2 r + L 2 -2R r Lcos(90° + 0 R \ 

R (R r +R,) 2 -L 2 
2(^Rj ■+■ R r + L sin 0 R ) 

R t =(R 2 r +L 2 + 2R R Lsin d R ) V2 . 

Substituting the range sum ( R T +R r ) in Eqs. 19 and 20 using Eq. 17 yields 


(18) 

(19) 

( 20 ) 
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( 21 ) 


r r = 


LQ-g 2 ) 

2e(\ + esmO R ) 


R t = 


L(c +1 + 2c sin 0 R ) 
2^(1 + ^ 8016 ^) 


( 22 ) 


Using the law of sines on the bistatic triangle defines the relationship between the range 
and angle values as 


R„ 


Rt 


sin(90 -6, ) sin(90 + 0 R ) sin /? 
Hence, the bistatic angle p is expressed as 


(23) 


P = sin 


P - sin 


( t _ n \ 



(24) 


(25) 


The direct and indirect method of measuring range sum (R T +R r ) applies to all 

target locations except when the target is in a forward-scattering configuration such that 
the target lies on the baseline between the transmitter and receiver. In the forward¬ 
scattering PBR configuration, (R t +R r ) = L and 0 R = -90 , making R r in Eq. 19 
indeterminate [2], 


In view of this limitation for PBR in forward-scattering configurations and 
inaccuracies arising from estimating the target’s location using a single PBR with 
erroneous time delay and AOA measurements, as well as inaccurate transmitter and 
receiver position information, the following sub-section introduces the least-squares 
solution for fusing bistatic triangle parameters from all PBR pairs. 


2. Least-Squares Intersection of Lines 

In a realistic scenario, AOA and range information for each PBR pair derived 
from measurements and solving trigonometric equations do not result in bearings 
intersecting at a single point (Figure 8). The least-squares solution derived in [32] finds a 
point that minimizes the sum of perpendicular distances from this point to all the lines. 
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This method (i.e., the one that solves the bistatic triangle parameters and estimates the 
target location by least squares bearing intersection) is referred to as the elliptical target 
location method. 



Legend: 

X Target 

- Txl bearing from Txl position and 0 T - Rx bearing from R r and 0 R for Txl 

- Tx2 bearing from Tx2 position and 0 T - Rx bearing from R r and 0 R for Tx2 


Tx3 bearing from Tx3 position and 0 T 


Rx bearing from R r and 0 R for Tx3 


Figure 8. Least-squares intersection of lines solution to three PBR case. 


A two-dimensional line is described by a point on the line 


a = 



(26) 


and its corresponding unit direction vector 


b = 



b r b = 1. 


(27) 
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The squared perpendicular distance from a point p to a line as illustrated in Figure 9 is 
expressed as 


D(p;a,b) = |(a-p)-((a-pfb)b| 2 
= (a - p) T (I - bb T )(a - p), 
and the sum of squared distance for K lines is 


(28) 


K 

D(p;a,b) = ^D(p;a.,b.) 

7=1 

K 


= E( a ;-P) r ( I - b ; b iK a i-P>- 


(29) 




Figure 9. Perpendicular distance from a point to a line. From [32], 


The corresponding objective function that finds the “best-fit” intersection point by 
minimizing the sum of squared distances for all lines is 

arg min 


P = 


£>(p;A,B). 


(30) 


Taking the derivative of the cost function with respect to p 


dD 

dp 


X-2(I-b ; b')(a_,-p) = 0, 


7=1 


Sp = q, 
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(31) 

(32) 


where 


S = £<I-b,bp, q = f](I-b,b> r 

' 1 ' ' (33) 

Solving for p in Eq. 32 gives the “best-fit” point of intersection 

P = S-'q. (34) 

For each PBR pair, two lines are defined after solving for their bistatic triangle 
parameters: the transmitter-to-target bearing and the target-to-receiver bearing. The 
transmitter-to-target bearing is defined by the transmitter position and bistatic triangle 
parameter^ while the target-to-receiver bearing is defined by the point determined by 

( 0 K , R r ) and 0 R . In a forward-scattering configuration, where R r is indeterminate, the 
target-to-receiver bearing is defined by the receiver position and the bistatic triangle 
parameter 0 R . 


3. Hyperbolic Target Location 

Apart from the elliptical target location method covered in Section II.C.l and its 
extension in Section II.C.2, a target’s position can also be estimated using a hyperbolic 
location technique. To derive an iterative least-squares method of estimating the target’s 
location given time delay measurement from the direct method illustrated in Figure 6, 
consider a two-dimensional PBR pair in a multi-bistatic radar scenario (Figure 10). 
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Figure 10. Single PBR in multi-bistatic radar scenario. 

Given the coordinates of the target (x e , y e ) , receiver position (x rx ,y rx ) and 
transmitter positions (x txi , y tx i ), the noiseless time delay measurement in Eq. 12 can be 


rewritten as 

c &Ti R Tj + R r L 

h (X„ X, u , X„) = p.-x^y + iy.-y^,) 2 
+yj(x e -x rx ) 2 +(y e - y rx ) 2 
-^(x rx -x tx4 ) 2 +(y rx -y txa ) 2 

where 


i = 

transmitter number i = 1,2, • • •, K , 

K 

number of transmitters, 

R T,i = 

I th transmitter-to-target range, 

Rr = 

target-to-receiver range, 

L 

baseline range, 
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AT = time delay between i th transmitter’s direct and indirect signal, 

X e = target’s position vector, [x e y e f, 

X txJ = i lh transmitter’s position vector, y txi ~J , 

X rx = receiver position vector, [x rx y rx ] 7 . 

When noisy time delay measurements are used to estimate the target location, Eq. 
35 is written as 

_v, = cAT + n t 

='>,(x„x w ,x„)+„, (37) 

where 

y, = i th noisy time delay measurement, 

AT = i th noiseless time delay measurement, 

n t = j th time delay measurement error. 


Since the function in Eq. 36 is a non-linear function of the target, receiver, and 
transmitter positions, the function h (X t ,, X (A ( , X rf j will be linearized by a Taylor series 


expansion about an initial estimate of the target’s location (x e0 , y e0 ). By retaining the first 
order terms, Eq. 37 can be written as 

y, - h, (X„, X ., X„ ) = 4 Ax, +Ay, + n ,, 

&, dy, (38) 

where 




+ 


x 


eO 


y e o. 


(39) 

For K time delay measurements (K PBR pairs), Eq. 38 can be represented by a 
linear model 


Y = H X + N , 

(Xxl) (K\2) (2x1) (Axl) 


( 40 ) 
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where 


Y = 


y,-A,(X, 0 ,X„.„X„) ' 

y 2 ->h(X„X„ s ,X„) 


, H = 


dJ\ 

dl\ 

dx e 

dy e 



dx e 


dh K 

dh K 

dx e 

dy e 


, X = 


Ax e 

a y. 


(41) 


The least-squares solution X that minimizes the sum of squares of difference between the 
measurements and the estimated function is defined as 

X = (H r H) H T Y, 

where 


(42) 


X = 


AXe 

Ay, 


The estimated target location in the current iteration is therefore 


x,=x e 0 +x. 


(43) 


(44) 


such that the target location estimate in the current iteration is used as the initial estimate 
X (0 in the subsequent iteration. 

D. MEASUREMENT AND LOCATION ERRORS 

The theoretical root mean square (rms) error dM of a radar measurement M can 
be expressed as [33] 

kM kM 


dM = 


y]2E/N 0 4 sTn 


(45) 


where k is a constant whose value is of the order of one, E is the received signal energy, 
and N 0 is the noise power per unit bandwidth. 
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For time-delay (range) measurements, k depends on the shape of the frequency 
spectrum S(f ), and Mis the rise time of the pulse. For angle measurements, £ depends 
on the shape of the aperture illumination A(x ), and M is the beamwidth. 


1. Time Delay (Range) Measurements 

The theoretical rms error in time delay measurements dT R for a rectangular pulse 

with pulsewidth r and limited by bandwidth B is approximately [34] 

/ _ V /2 


dT„ 


k 4BE/N„ j 


and can be expressed as a function of SNR: 


dT D 


1 


b4sTn' 

This assumes Br »1, which is not always satisfied, 
rectangular pulse ( Bt ~ 1) for which [33] 

r 


dT„ 


2.1\fs 1 ~N 


(46) 


(47) 

A more accurate model is a quasi- 


(48) 


2. Angle Measurements 


The theoretical rms error of AOA measurements dO R for an antenna with uniform 


(rectangular) amplitude illumination across the aperture is [33] 


dd R =- 


x/31 


0.628# e 

V2 


ttD(2E/N 0 ) (2 E/N 0 ) 

and can be expressed as a function of SNR: 

°. 628 ^ 


0 B = 


4sTn ’ 
0 . 88/1 
D 


where D is the antenna dimension and 0 B is the half-power beamwidth. 


(49) 


(50) 

(51) 
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3. Transmitter and Receiver Position Accuracy 

Transmitter positions on merchant ships are updated by differential GPS (DGPS) 
systems and made available to surrounding vessels by onboard AIS units. There are two 
types of AIS transceivers [4]: 

1. Class A onboard commercial vessels 

2. Class B, used by leisure and smaller crafts 

For the purpose of this thesis, merchant ships are assumed to be equipped with Class A 
AIS that broadcasts the vessel’s unique identification, position, course, and speed 
information every 2 to 10 seconds while underway, and every three minutes while at 
anchor at a power level of 12.5 W [30], 

Receiver positions on warships are also provided by onboard DGPS systems. In 
this thesis, the U.S. Coast Guard’s (USCG’s) DGPS service accuracy of 2 distance-root- 
mean-square (drms) [35] is used to model transmitter and receiver position accuracies. 


4. Receiver-to-Target Range Error 

As explained in Section C.l, the receiver-to-target range R r is calculated from the 
range sum (Rj. + R r ), receiver look angle 6 R , and baseline L measurements. Assume 
that the measurement errors associated with range sum (R T +R r ) , receiver look angle 0 R , 

and baseline L are uncorrelated, zero-mean random processes having Gaussian 
distribution with standard deviation equal to measurement rms error. The geometry 
dependent root-sum-squared (rss) error of R r corresponding to Eq. 19 is expressed as [2] 


d R r = 


8R b 


-|2 


d(R T +R r ) 


d (R t +R r ) 


+ 


SR, 

dL 


d L 


+ 


SK 

yd& R 


d0 B 


(52) 


where d (R t +R r ) , dL and d 0 R are the rms errors for (R t +R r ) , L , and 0 R , 
respectively, and determined by Eq. 45. 

From Eq. 52, the components of the rss error estimate of R r can be expressed as a 
function of eccentricity e and receiver look angle 0 R [2]: 
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SR r 

1 + 6 + 2e sin 0 R 


8(R t +R r ) 

2(1+ e sin 0 R ) 2 

(53) 

8R r {c +1) sin 0 R + 2e 


dL 

2(l + esin 0 R ) 2 

(54) 

$ 

1 

1 

L(l-e 2 )cosO R 


de R 

2(1+ esin 6 R f ’ 

(55) 


where the elliptical iso-range contour eccentricity is e = L/(R T +R R ) and each partial 

derivative component defines the slope of the error surface with respect to each 
measurement variable. 

E. ERROR ELLIPSE PARAMETERS 

The error ellipse provides a graphical means of viewing uncertainty associated 
with position estimates. The error ellipse is described by three parameters: (1) semi-major 
axis, a ; (2) semi-minor axis, b ; and (3) orientation, 0 (Figure 11). 


v 
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The error ellipse parameters are described by the covariance matrix of the target’s 
position estimate, C x , given by 


C „ = 


cr p a cr 

x e ^ x e y e x e y e 

p <j <j g 2 

^y e x e y e x e y e 


where cr" , g 2 and p v v are the variance and correlation associated with the target’s 

'X'e ye y e^e 

location estimate. 

For the hyperbolic target location method, the target location estimate covariance 
matrix is expressed as 

C x e y e = J vv e ] ’ ( 57) 

where C„ is the covariance matrix for time delay measurement errors and J, is the 

n J x e T e 

Jacobian matrix of partial derivations in Eq. 57: 

X 0 0 0 ' 

r _ 0 0 0 

" 0 o 0 


) 0 


dh t 

d\ 

dx e 

dy e 

dh 2 

dh 2 

dx e 

to. 

dh K 

dh K 

dx e 



The semi-major axis a and semi-minor axis b are expressed as a function of 


elements in C „ as [36] 

x e.ye 


a = \h( a ^ + ) + \U K + ) 2 + ( p y^ a y, a 7 ) 2 ’ 


b = A -(a x +CTy)-. -(a x +a x ) + (p cr cr ) . 

\ O V ye / \i A \ x e / \/ Te x e / 
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which can be further reduced to a general form 


a = 
b = 


'JzLA 

sfzlA 


/l, >^2 


(62) 


where * s the critical value of the chi-square distribution with 2 degrees of freedom 


and evaluated at a (Figure 12). 
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Figure 12. Chi-square probability density function with 2 degrees of freedom. 
The area to the right of x 1 critical value is a . 


There are two forms of equations used to compute the orientation of the error 


ellipse [36, 37]: 


1 


0 = — tan -11 
2 


^ 2p „ cr, (j ^ 




6 = — tan 1 
2 


<r x . 

2 


ye J 

7 2 A 


2 p, 


<7 <7 

x e y e x e y e j 


(63) 

(64) 


30 



where Eq. 63 and Eq. 64 represent clockwise and anticlockwise rotation of the error 
ellipse to achieve statistical independence (Figure 13). 



Figure 13. Error ellipse rotation to achieve statistical independence. From [38]. 


To understand the relationship of Monte Carlo simulations of target locations to 
their corresponding error ellipses, see the derivation of Eq. 62 and Eq. 63 in Appendix A. 
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III. MATLAB AND FEKO MODELING 


A. PROBLEM SETUP 

The maritime scenario with a single receiver and multiple transmitters depicted in 
Figure 2 will be employed in all models. The problem setup assumes a warship (receiver) 
equipped with broadband EW receiver and DF capability deployed to monitor maritime 
traffic flow at the Straits of Singapore. The targets are assumed to be low-RCS targets 
with infrequent transmissions and/or operating with LPI radars while civil marine radars 
(CMRs) on merchant ships will serve as opportunistic illuminators in PBR configuration 
[4], Regulation 19 of SOLAS Chapter V requires that all merchant ships and warships 
carry automatic identification systems (AISs) which share information on the ships’ 
identity, position, course, speed, navigation status, and safety-related information [30]. 

The target detection process starts with the warship receiving AIS transmissions 
from surrounding transmitters, allowing the warship to form a tactical picture with 
knowledge of the transmitter positions and the corresponding baseline ranges. Onboard 
EW and DF receivers will identify merchant ships, which will allow separation of direct 
and indirect (target echo) transmissions from the same CMR. The time delay between 
reception of the direct and indirect transmissions, DF of target echo and baseline range is 
used to compute the target-to-receiver range in all PBR configurations except in a 
forward-scattering configuration. An extension to this elliptical target location method is 
to define transmitter-to-target and target-to-receiver bearings for each PBR and estimate 
the target location by finding a “best-fit” intersection via least-squares. An alternative to 
the abovementioned target location method is the hyperbolic target location approach, 
which uses time delay measurements and knowledge of the transmitters and receiver 
position to estimate the target’s position iteratively. 

In this thesis, technical specifications from commercially available systems are 
used to model transmitters, EW and DF receivers. Technical specifications of UK-based 
Kelvin Hughes’ MantaDigital Radar (see Appendix B) in the S-band and X-band are used 
to model the opportunistic transmitters. 
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For the purpose of this thesis, merchant ships are assumed to be equipped with 
Class A AIS that broadcasts the vessel’s unique identification, position, course and speed 
information every 2 to 10 seconds while underway, and every 3 minutes while at anchor 
at a power level of 12.5 W [30]. Receiver positions on warships are also provided by 
onboard DGPS systems. For the simulation, the USCG’s DGPS service accuracy of 2 
drms [35] is used to model transmitter and receiver position accuracies. 

Onboard the warship, the EW receiver would measure and extract transmission 
parameters, after which emitter classification and identification takes place. Integrated 
with the EW receiver is the DF receiver which determines the AOA of transmission 
signals. The EW receiver specifications used in the models are based on UK-based 
Teledyne Defence Limited’s high performance QR026 receiver (see Appendix C). This 
wideband EW receiver operates from 0.5 GHz to 18 GHz and has a sensitivity of -70 
dBm. The DF receiver specifications are based on South Africa’s Poynting Innovations 
(Pty) Ltd’s DF-AA0062 antenna (see Appendix D), which has coverage from 20 MHz to 
6 GHz. 


B. FEKO MODEL 

The EM simulation software FEKO was employed to model and analyze the RCS 
scattering properties of the maritime target. A generic frigate is chosen as the maritime 
target. A three-dimensional model of the frigate was created in Computer Simulation 
Technology (CST) Microwave Studio (MWS) in [4] and imported into FEKO where free 
space monostatic and bistatic RCS returns at S-Band and X-Band can be generated. 
Figure 14 shows the FEKO mesh model and a picture of the actual frigate on which the 
model is based, while Figure 15 shows the model’s coordinate system. All surfaces are 
assumed to be perfect conductors. 
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Figure 14. Frigate FEKO model (top) and actual RSN Formidable-class frigate 

(bottom; from [39]). Side-profile. 

Figure 16 shows the frigate monostatic RCS and bistatic RCS at 10° incident 
angle at 3.05 GHz. The difference in monostatic and bistatic RCS can vary up to 70 dB 
for different look angles. The frigate’s bistatic RCS returns vary with bistatic angle (3 
and is larger than the monostatic RCS in the forward scatter direction. The PBR 
configuration capitalizes on this characteristic to detect low-RCS targets using 
opportunistic transmissions. 
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Figure 15. Frigate FEKO model and coordinate system. 
Frigate Monostatic RCS Frigate Bistatic RCS 



Figure 16. Frigate’s monostatic RCS (dBsm) (left) and bistatic RCS (dBsm) 
with incident angle of 10° (right) at 3.05 GHz. 

Figure 17 shows the frigate’s bistatic RCS returns at 1° and 0.1° resolution, 

respectively, against the azimuth return angle. From the plots, it is evident that small 

variations in azimuth return angles may result in significant difference in RCS returns. 

Simulations using the same model in FEKO also suggest that small variations in the 

azimuth incident angle will result in significant variations in RCS returns. Therefore, to 
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ensure that accuracy of the model is maintained, it would be ideal to call FEKO to 
compute the exact RCS returns given the required incident and receive angles. 


Frigate Bistatic RCS 



150 iao 210 
Phi [deg] 


Figure 17. Frigate bistatic RCS (dBsm) with incident angle of 10° at 3.05 GHz. 
RCS at 1° resolution (left) and 0.1° resolution (right). 


C. MATLAB MODEL 

The maritime scenario is modeled in MATLAB with the frigate target RCS from 
FEKO. Detection coverage and target location methods are simulated using the 
MATLAB models, after which the results are examined to propose optimal PBR 
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configurations and geometry, as well as to investigate their effects on estimated target 
location uncertainty. 

1. Detection Coverage Model 1 

The original MATLAB detection coverage model (.Detection Coverage Model 1 ) 
computes the received SNR given the target’s position information. This model requires 
the following input file and parameters to generate a detection coverage summary plot: 

• Target path information in an Excelfile 

• Target CADFEKO model file (.cfx file) 

• Target default EDITFEKO file (.pre file) 

• Transmitter parameters 

• Receiver parameters 

The target’s path information Excelsheet is of the form given in Table 3, where 
the model is required to execute a FEKO run for each PBR pair at each time step. The 
target’s position and velocity vector at each time step can represent the actual target’s 
course or it could be a grid of equally spaced arbitrary target positions. The detection 
coverage model generates the received SNR along the target’s path for the former, and a 
detection coverage plot (receiver SNR are grid points) for the latter. 


Table 3. Target path information format in Excelsheet. 


Time 

Position Vector 


Velocity Vector 


h 

x 1 

yi 

Zi 

v ,i 

V yl 
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V y2 

V z2 
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X n 

y* 

Zn 

V xn 


V z« 


The target’s CADFEKO and EDITFEKO files are required for FEKO to model 
and generate the RCS returns from the target for given incident and receive angles. The 
CADFEKO file contains the target’s mesh model while the EDITFEKO file includes 
execution parameters such as frequency, incident angle, and receiver angles. For each 
FEKO run, parameters in the EDITFEKO file are modified prior to running the FEKO 
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Solver function on the CADFEKO model. The RCS returns are then read from the .out 
file. The transmitter, EW, and DF receiver S-band and X-band parameters are based on 
technical specifications of commercially available systems. 

Figure 18 gives an overview of Detection Coverage Model 1 for each transmitter- 
target-receiver pair at each time step. The entire MATLAB model and its interface with 
Exceland FEKO are automated. The model begins by computing the incident and receive 
angles at the target given the transmitter’s, receiver’s, and target’s current positions. 
Given the incident and receive angles, the EDITFEKO file is modified to indicate the 
required signal incidence and receive angles. The FEKO solver is then executed on the 
target’s CADFEKO model to generate an .out file containing the RCS return for the 
current transmitter-target-receiver geometry. 



Figure 18. MATLAB Detection Coverage Model 1 flowchart. 


2. Detection Coverage Mlodel 2 

Figure 19 gives an overview of the Detection Coverage Model 2 for each 
transmitter-target-receiver pair at each time step. Unlike Model 1, Model 2 extracts RCS 

returns via interpolation given bistatic incident and receive angle from a pre-computed 
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RCS table of the form given in Table 4. This reduces the simulation duration greatly as 
the model does not need to execute FEKO for each angle. 



Figure 19. MATLAB Detection Coverage Model 2 flowchart 


The Detection Coverage Model 2 also computes the received SNR given the 
target’s position. This model requires the following input file and parameters to generate 
a detection coverage summary plot: 

• Target path information in an Excelfile 

• Target RCS table in an Excelfile 

• Transmitter parameters 

• Receiver parameters 
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Table 4. RCS table format in Excelsheet. 



The target’s path information Excelsheet and RCS table would be of the form 
given in Table 3 and Table 4, respectively. Like Model 1, the target’s position and 
velocity vector at each time step can either represent the actual target’s course or it could 
be a grid of equally spaced arbitrary target positions. The detection coverage model 
generates the receive SNR along the target’s path for the former, and a detection 
coverage plot (receiver SNR are grid points) for the latter. The table of RCS values for 
various incident and receiver bistatic angles is pre-computed in FEKO using the target’s 
CADFEKO model. The bistatic angle resolution for the RCS table can be improved by 
generating a higher resolution RCS table prior to simulating detection coverage. Lastly, 
transmitter, EW, and DF receiver S-band and X-band parameters are based on technical 
specifications of commercially available systems. 

3. Simulation Duration and Accuracy Trade-off 

Preliminary simulations using Detection Coverage Model 1 presented the issue of 
simulation duration and accuracy trade-off. Referring to the flowchart of Model 1 in 
Figure 18, it takes approximately 4 seconds for FEKO to compute the RCS each time the 
FEKO Solver is called. The time taken is acceptable if the model is required to execute 
FEKO via the model only a few times. However, in generating detection coverage plots, 
the frigate path Excelfile contains equally spaced grid points across an area of interest. In 
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this case, the model calls FEKO at each grid point for each PBR pair. For example, 
consider a 40 km x 40 km area with grid points spaced 1 km apart (Figure 20); there will 
be 41 x 41 = 1681 grid points. If there is one opportunistic transmitter in the grid area, it 
will take approximately 1681 x 4 x 1= 6728 seconds =1.9 hours to generate the detection 
coverage for the grid area. However, if there are 8 opportunistic transmitters in the grid 
area, it now takes approximately 1681 x 4 x 8 = 53,792 seconds = 14.9 hours to generate 
the detection coverage. 
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Figure 20. Detection coverage of a 40 km x 40 km area of interest grid points at 

1 km resolution. 


The large number of simulation cases, compounded with the need to simulate 
various grid sizes/resolutions and number of transmitters, resulted in the alternative 
detection coverage model (Detection Coverage Model 2). 

Referring back to Figure 17, small variations in return angles may result in 
significant differences in RCS return; therefore, Model 2 reduces detection coverage 
simulation duration at the expense of accuracy due to low-RCS resolution. 
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Figure 21 illustrates the simulation duration and accuracy trade space for Model 1 
and Model 2. Assume that the RCS table used by Model 2 has bistatic angle resolution of 
1 degree and its current position in the simulation trade space is indicated in Figure 21. 
The accuracy of Model 2 can be improved with minimal increase in simulation duration 
by pre-computing RCS tables with higher resolution bistatic angles. Generating a high 
resolution RCS table requires longer computing duration compared to a lower resolution 
RCS table. However, the additional computing duration occurs one time upfront and 
greatly reduces the simulation duration of Detection Coverage Model 2. 

In view of the constraint on computing speed, Detection Coverage Model 1 and 
Detection Coverage Model 2 will be generated but Model 2 will be used to simulate 
target detection coverage and target location performance for the rest of this thesis. 



Figure 21. Simulation duration and accuracy trade space for detection coverage 
models. 
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4. Target Location Model 

The MATLAB target location model estimates the target’s location given a set of 
erroneous measurements and parameters. This model requires the following input file and 
parameters: 

• Target path information in an Excelfile 

• Erroneous time delay, DF measurements 

• Erroneous transmitter and receiver positions 

The target’s path information Excelsheet would be of the form given in Table 3. 
The target’s position and orientation are used to generate measurements and compute the 
SNR at the receiver, which in turn is used to model rms errors in the measurements. The 
theory behind modelling measurement errors is covered in Section II.D. For the 
hyperbolic target location method, the error ellipse associated with each position estimate 
is computed. 

D. MODEL VERIFICATION 

Prior to using the Detection Coverage Models to generate simulation results, it is 
necessary to verify the model against bistatic radar principles and/or using simple targets. 
Here, the detection coverage models were verified using the following approaches: 

• Target with constant RCS 

• Target plane at various orientations 

According to Eq. 2 and Eq. 3 the SNR of the target echo depends on transmitter 
parameters, receiver parameters, target RCS o B , transmitter-to-target range R T , and 

target-to-receiver range R r . When transmitter parameters, receiver parameters, and target 
RCS are kept constant, the SNR depends solely on R T and R r . For each pair of constant 
R t and R r ranges, the contour of constant SNR traces out an oval of Cassini, as covered 
in Section II.B. 

To verify that the Detection Coverage Models satisfy bistatic radar principles, 
consider a target whose RCS is constant regardless of aspect angle or range. The resultant 
detection coverage should exhibit elliptical SNR contours, as depicted in Figure 5. The 
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detection coverage and contour plot generated from the Detection Coverage Models is 
shown in Figure 22. The results agree with bistatic radar theory such that the SNR 
contour values satisfy Eq.2. 


Detection Coverage SNR Contours (dB) 



Figure 22. Detection coverage contours for constant RCS target. 

Apart from verifying the model against bistatic radar principles, it was tested 
using a simple plate target at various orientations. The FEKO model of the target plate 
used is shown in Figure 23. Figure 24 shows the S-band RCS for the target plane where 
the angle represents the angle (f) from the x-axis horizontal in the xy-plane (6 = 90 ) of 
the FEKO model. 
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Figure 23. Target plate FEKO model. 
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Figure 24. Target plate S-Band azimuth RCS (dBsm) at cp=90° incidence. 
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To verify that the MATLAB Detection Coverage Model is sound for various 
target orientations, detection coverage plots for the target plate are generated at several 
rotations of the plate around the z-axis. From the plate’s RCS in Figure 24, one would 
expect higher RCS and hence SNR at the faces of the plate, that is, at 90° and 270°. In 
addition, as the plate’s RCS returns exhibit symmetry about a certain axis, the detection 
coverage plots at selected target orientations demonstrate symmetrical characteristics. See 
Appendix E for the target plate detection coverage plots and discussion of its RCS. 

The elliptical and hyperbolic target location methods outlined in Section II.C were 
verified using measurements in the absence of errors. Using error-free measurements, 
estimated target locations have zero errors for both target location methods. 

The hyperbolic target location’s associated error ellipse semi-major length, semi¬ 
minor length, and orientation are verified against Monte Carlo simulations. For example, 
it was verified that approximately 90% of the Monte Carlo estimates fall within a 90% 
confidence ellipse. 
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IV. SIMULATION RESULTS 


In this chapter, the detection coverage for various PBR transmitter-target-receiver 
geometries is examined using Detection Coverage Model 2, and PBR configurations and 
geometry for optimal detection coverage are proposed. The results for elliptical and 
hyperbolic target location methods and their corresponding uncertainty ellipses are 
presented. This chapter begins by investigating how the number of opportunistic 
transmitters, their ranges, and geometry relative to the target and receiver affects 
detection coverage in an area of interest. Detection results from Detection Coverage 
Model 1 and Model 2 along the target’s arbitrary path are also presented. Target location 
results from elliptical and hyperbolic target location methods are then presented with 
their corresponding error ellipse. The effects of transmitter-target-receiver geometry on 
each target location method’s spatial error distribution are explored. All simulation plots 
and results are included in Appendices F through H. 

A. DETECTION COVERAGE 

1. Bistatic RCS 

To better understand the detection coverage plots of the PBR scenario presented 
in this section, consider the bistatic RCS of the target (Figure 14) with incident angle 10° 
at 3.05 GHz and 9.41 GHz in Figure 25. The frigate’s bistatic RCS returns vary with the 
incident and receive angles with peaks occurring at certain return angles. Overall, they 
have similar return patterns in S-band and X-band. It is also observed that RCS will have 
a peak in the forward-scattering look angle. 

Detection coverage plot results in this section are obtained using Detection 
Coverage Model 2 with a RCS table of 1° angle resolution. The RCS values from the pre¬ 
computed table are extracted by interpolation if incident and/or receive angles are not 
integers. Each point (x, y) on the detection coverage plot represents the total received 

SNR at the receiver when the target is at position (x, y). The colour at each point 
indicates the value of total SNR at the receiver. 
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Frigate Bistatic RCS at 3.05 GHz 

90 


Frigate Bistatic RCS at 9.41 GHz 

90 


Figure 25. Target bistatic RCS (dBsm) with incident angle 10° at 3.0 GHz (left) 

and 9.41 GHz. 


2. General Observations 

Referring to the detection coverage results (see Appendix F), the SNR of a target 
echo at the receiver strongly depends on the bistatic RCS characteristic of the target. 
Higher SNR will always be observed when the target-to-receiver LOS lays in the peak of 
the target’s bistatic RCS returns. Since the target’s bistatic RCS exhibits high RCS 
returns in the forward-scattering look angle (Figure 25), a PBR with the transmitter- 
target-receiver in the forward-scattering configuration results in the best detection at the 
receiver for a fixed range. It is also evident from the detection coverage results that 
transmitter-target-receiver in the back-scattering configuration would result in good 
detection with target echo SNR decreasing as the target’s range to receiver increases. 
Figure 26 highlights target positions in the forward and back-scattering configuration on 
a detection coverage plot where the transmitter is positioned at (-5000,0) and receiver 
at(0,0). 
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Figure 26. Target positions in forward-scattering and back-scattering 
configuration on detection coverage plot. 


From the detection coverage results (see Appendix E), it was observed that target 
echo SNR is generally high when the target is within the “net” of transmitters-receiver 
configuration with the highest SNR occurring when the PBR is in the forward-scattering 
configuration. This observation agrees with the findings in [4] which state that the 
optimal transmitters-receiver constellation is to have transmitters at every quadrant with 
respect to the receiver. Outside the “net” of transmitters-receiver, detection is usually 
better when the PBR is in the back-scattering configuration. Across all PBR geometries, 
detection gaps—areas with the low target echo SNR at the receiver—occur when the 
target is outside the transmitter-receiver cluster or away from bearings extending from the 
transmitter-receiver baselines (Figure 27). 
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Detection Coverage SNR (dB) 



Figure 27. Detection gaps and bearings extending from Tx-Rx baseline in 

detection coverage plot. 


3. Vary Number of Transmitters 

The number of available opportunistic transmitters affects the detection coverage 
area and strength of the target echo at the receiver. Figure 28 shows that overall detection 
coverage improves as the number of opportunistic transmitters increases. The result 
presented in this thesis assumes that target echo is received from all transmitters shown in 
the detection coverage plot. The percentage of grid points with received SNR greater than 
10 dB within a 10 km radius from the receiver for a different number of transmitters 
around the receiver is given in Table 5 and 0. For a target orientation of 0°, detection 
coverage increases from 17.0% to 61.6% as the number of opportunistic transmitters 
increase from one to eight (Table 5). 
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Figure 28. Detection coverage plots at S-Band for with target plate at 0° 
orientation as the number of transmitters varies. 

The detection statistics presented in Table 5 and 0 suggest that S-band 
opportunistic transmitters result in better detection than X-band opportunistic transmitters. 
This is attributed to the variation in RCS returns in X-band compared to S-band. 
Referring to Figure 25, peak RCS at X-Band is higher, but narrow lobes result in low 
probability of intercepting a lobe. Nevertheless, detection coverage patterns are consistent 
between S-band and X-band transmissions. As the difference in detection coverage is not 
significant in detection coverage plots, subsequent sections of this thesis are focused on 
results from S-band transmissions. 
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Table 5. Percentage of grid points with SNR greater than 10 dB within a 10 
km radius centered at the receiver. Target orientation at 0°. 


Target at 0° orientation 

# of Tx 

% SNR > 10 dB in 10 km radius 

S-Band 

X-Band 

1 

17.0% 

15.1% 

2 

33.8% 

29.5% 

3 

42.0% 

36.4% 

4 

47.5% 

40.3% 

8 

61.6% 

54.8% 


Table 6. Percentage of grid points with SNR greater than 10 dB within a 10 
km radius centered at the receiver. Target orientation at -90°. 


Target at -90° orientation 

# of Tx 

% SNR > 10 dB in 10 km radius 

S-Band 

X-Band 

1 

16.4% 

15.7% 

2 

33.4% 

29.5% 

3 

40.0% 

35.4% 

4 

47.5% 

40.3% 

8 

61.6% 

54.8% 


4. Vary Transmitter Range 

The transmitter to receiver range determines the size of the transmitter-receiver 
“net” and forward-scattering area. Detection coverage plots in Figure 29 suggest that the 
area of detection coverage increases as the transmitter to receiver ranges increase; 
however, the highest target echo SNR in the grid decreases with increase in range. This 
observation is consistent with the findings in [40] with regards to sensitivity dependency 
on transmitter range and transmitter-receiver geometry. The percentage of grid points 
with received SNR greater than 10 dB within a 20 km radius from the receiver- 
transmitter range of 5 km, 10 km, and 15 km are listed in Table 7. For target orientation 
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of 0°, detection coverage increases as the transmitter to receiver range increases from 5 
km to 10 km. However, the detection coverage decreases when the range further 
increases to 15 km unless more than four transmitters are present. 



Figure 29. Detection coverage plots at S-Band with target plate at 0° orientation 

as the transmitter range to receiver increases. 
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Table 7. Percentage of grid points with SNR greater than 10 dB within a 20 
km radius centered at the receiver. Target orientation at 0°. 



% SNR > 10 dB in 20 km radius 

Tx-Rx Range 

# of Tx 

5 km 

10 km 

15 km 

1 

7.3% 

8.7% 

5.9% 

2 

14.4% 

17.3% 

12.7% 

4 

19.7% 

28.8% 

19.9% 

8 

31.2% 

37.1% 

38.0% 


5. Vary Transmitter-Target-Receiver Geometry 

Two types of transmitter-target-receiver geometries were explored: (1) receiver- 
centered geometry, as seen in Figure 30, and (2) transmitter-clustered geometry, as seen 
in Figure 31. In a receiver-centered geometry, all observations made earlier in this 
chapter hold. High target SNR is received when the target is in the transmitter-receiver 
“net” with the best detection occurring when the transmitter-target-receiver is in the 
forward-scattering configuration. Outside the transmitter-receiver “net,” detection is the 
best when the transmitter-target-receiver is in the back-scattering configuration. Apart 
from areas that satisfy a back-scattering configuration, it is observed that there are 
additional areas that exhibit good detection coverage as a result of the target’s bistatic 
RCS patterns. Multiple PBR capitalizes on this to allow detection of low-RCS targets. 

In a transmitter-clustered geometry, good detection is experienced when the target 
lies in the cluster of transmitters or in the beam extending from the receiver to the cluster 
of transmitters. Outside the transmitter-receiver cluster, detection is slightly better at 
areas that satisfy back-scattering configuration and “lines” that are shaped like a 
hyperbolic curve—hyperbolic contours are typical of time difference of arrival (TDOA) 
measurements [37]. However, for a four opportunistic transmitter’s scenario, the received 
target’s echo is still below 10 dB outside the transmitter-receiver cluster. 
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Figure 30. Detection coverage at S-Band for 4 Tx in receiver-centered geometry 

and target plate at 30° orientation. 
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Figure 31. Detection coverage at S-Band for 4 Tx in transmitter-clustered 
geometry and target plate at 30° orientation. 
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Up to this point, all transmitters in a scenario were placed in a predictable fashion 
within an area of interest. The next series of detection coverage plots randomly positions 
20 transmitters at a range of 5 to 20 km from the receiver. Earlier observations in regard 
to areas of good detection coverage are evident in the detection coverage plots. However, 
strips of detection gaps—positions with the lowest target echo SNR at the receiver—are 
evident in Figure 32, Figure 33, and Figure 34. 



Distance (m) x 

Figure 32. Detection coverage at S-Band for 20 Tx randomly positioned around 
Rx at 5-30 km range and target plate at 0° orientation. 

The direction and location of the detection gaps are determined by the target’s 
velocity vector and hence its orientation. Referring to Figure 32 where the target velocity 
vector is horizontal and points to the right and, hence, the target orientation is 0° across 
all grid points, detection gaps will present themselves as horizontal strips. Applying this 
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explanation for target orientations of -90° and 45°, detection strips would be vertical and 
angled at 45°, respectively (Figure 33 and Figure 34). The strips of detection gaps are 
observed to be present when the target’s heading is towards the transmitters-receiver 
cluster. 


Detection Coverage SNR (d8> 



Distance (m) 


Figure 33. Detection coverage at S-Band for 20 Tx randomly positioned around 
Rx at 5-30 km range and target plate at -90° orientation. 
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Figure 34. Detection coverage at S-Band for 20 Tx randomly positioned around 
Rx at 5-30 km range and target plate at 45° orientation. 

6. Target Path Detection 

Detection Coverage Model 1 and Model 2 compute detection performance given 
the target’s path, with Model 1 using more accurate RCS values from FEKO and Model 2 
extracting RCS values from a pre-computed table given bistatic incident and receive 
angles. The trade-off between accuracy and simulation duration illustrated in Section 
III.C.3 limits the use of Model 1 as a tool to study detection performance. Nevertheless, 
target path detection performance results from both models are included in Appendix G. 

Five arbitrary target paths in Figure 35 are used to generate target path detection 
performance as the target moves through the area of interest. The received target echo 
SNR as the target transverses through the area of interest is given in Figure 36. The 
detection performance results here align with previous observations, where detection 
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coverage is good within the transmitters-receiver cluster and in the back-scattering 
configuration outside the transmitter-receiver cluster. 
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Figure 35. Five target paths used to generate target path detection performance 

results. 
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Figure 36. Detection coverage along target’s path at S-Band for 4 Tx (left) and 
20 Tx (right). RCS extracted from pre-computed RCS table. 
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B. TARGET LOCATION ESTIMATION 


1. General Observations 

To examine target location estimation accuracy for elliptical and hyperbolic target 
location methods, results from both methods are presented in target location error color 
plots and estimated target location scatter plots. In a target location error plot, the color at 
each point (x, y ) represents the error in location estimate when the true target is located 

at (x, y). For estimated target location scatter plots, the estimated target location for each 

round of Monte Carlo simulation is marked such that spatial error distribution for each 
location method can be investigated. 

Referring to the target location error color plots where measurement errors are 
independent of SNR (see Appendix H), the hyperbolic location method generally 
estimates target location more accurately compared to the elliptical location method as it 
uses fewer measurement sources. The hyperbolic method uses time delay measurements 
and transmitters and receiver position data, while the elliptic method uses DF 
measurements, time delay measurements and transmitter and receiver position data. 
When measurements with errors independent of SNR are used, it was observed that the 
effects of dilution of precision on target location estimates are less pronounced for the 
hyperbolic method (Figure 37). 

To better model errors present in measurements used to estimate the target’s 
location, target location models will include SNR-dependent errors in the DF and time 
delay measurements, as outlined in Section II.D. Simulating target location accuracy with 
various combination of erroneous and error-free measurements suggest that DF 
measurement is the dominant source of error in the elliptical method, while time delay 
measurement is the dominant source of error in the hyperbolic method. It was also 
observed that the effect of 2 meters drms error in transmitter and receiver position data 
are negligible compared to SNR-dependent DF and time delay measurements. 
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Figure 37. Dilution of precision comparison between elliptical and hyperbolic 
methods using SNR-independent measurement errors. 


Referring to the target location error results where measurements with SNR- 
dependent measurement errors are used (Appendix H), the hyperbolic method, which 
uses fewer sources of measurement and data, produces more accurate location estimates 
compared to the elliptical method. However, overall location accuracy is poor due to low 
SNR and large DF and time delay measurement errors. 

The results from the target location accuracy color plots mirror the observation 
made in detection coverage plots. Target location accuracy is best when the target lies 
within a transmitters-receiver cluster in receiver-centered geometry (Figure 38). 
Increasing the transmitter to receiver range effectively increases the transmitters-receiver 
cluster area but reduces the overall location accuracy within the cluster. 
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Target Position Estimate Error (m) 



Distance (m) x iq 4 

Figure 38. Target location errors using hyperbolic method (S-Band) with SNR- 

dependent measurement errors. 

2. Error Ellipse of Target Position Estimate 

To examine the spatial error distribution of target position estimates, 100 Monte 
Carlo simulations are performed for each transmitters-target-receiver geometry. The 
target position estimate from each Monte Carlo simulation is marked on a two- 
dimensional plot to produce a target location scatter plot. For the hyperbolic target 
location method, parameters of each estimate’s error ellipse can be calculated for various 
levels of confidence. The equations used to compute the error ellipse semi-major axis 
length, semi-minor axis length, and orientation is covered briefly in Section II.E and 
derived in Appendix A. 

Transmitter-receiver geometries can be broadly classified into two groups: (1) 
target within Tx-Rx cluster and (2) target outside Tx-Rx cluster. 
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Figure 39 and Figure 40 show the target location scatter plot from the elliptical 
and hyperbolic methods respectively when group (1) geometry occurs. The estimates 
from all 100 Monte Carlo simulations cover an approximately circular area, suggesting 
that the extent of target estimate uncertainty in the x and y direction is similar. The area 
covered by the scatter points is referred to as the error ellipse for position estimate. The 
error ellipse of estimates from the hyperbolic method is smaller than that from the 
elliptical method. This suggests that the hyperbolic method estimates target position more 
accurately. 

Similarly, Figure 41 and Figure 42 show the target location scatter plot from the 
elliptical and hyperbolic methods respectively when group (2) geometry is considered. 
The estimates from all 100 Monte Carlo simulations in this case form an ellipse, 
exhibiting DOP. The resultant error ellipse from using the elliptical and hyperbolic 
method differs greatly with their semi-major axis appearing approximately orthogonal. 
To understand how measurement errors translate to uncertainty in the target location 
estimates, consider the graphical illustration in Figure 44 through Figure 47. 


Estimated Target Locations using Elliptical Method 


X Tgt 
Tx 
Rx 

• Estimates 


5000 r 

4000 

3000 

2000 

— 1000 
CD 

c 0[ 

£ 

GO 

b -1000 
-2000 
-3000 
-4000 
- 5000 ' 

-5000 -4000 -3000 -2000 -1000 0 1000 2000 3000 4000 5000 

Distance(m) 

Figure 39. Target position estimate scatter plot from elliptical method for target 

in Tx-Rx cluster. 
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Estimated Target Locations using Hyperbolic Method 
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Figure 40. Target position estimate scatter plot from hyperbolic method for 

target in Tx-Rx cluster. 
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Figure 41. Target position estimate scatter plot from elliptical method for target 

outside Tx-Rx cluster. 
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Estimated Target Locations using Hyperbolic Method 
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Figure 42. Target position estimate scatter plot from hyperbolic method for 

target outside Tx-Rx cluster. 

In the elliptical location method, bistatic triangle parameters are computed using 
DF and time delay measurements as well as transmitters and receiver position data. The 
transmitter-to-target range R, , receiver to target range R r , and DF measurements are 

used to define several bearings which would intersect at the target’s true position if all 
measurements are error-free. However, in reality, these bearings do not intersect at a 
single point, hence the elliptical method finds “best-fit” point of intersection and returns 
this point as the estimated target position. Graphically, uncertainty associated with each 
transmitter-to-target or receiver-to-target bearing is defined as a sector extending from 
the transmitter or receiver towards the target’s position. The intersection of all sectors 
approximates the uncertainty associated with the target position. For geometry (1), the 
target position estimates from the elliptical method results in a relatively circular and 
small error ellipse (Figure 43), while for geometry (2), the error ellipse is eccentric with 
the semi-major axis aligned with the receiver’s LOS (Figure 44). 
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Figure 43. Uncertainty area associated with position estimates from elliptical 

method for target in Tx-Rx cluster. 



Figure 44. Uncertainty area associated with position estimates from elliptical 

method for target outside Tx-Rx cluster. 
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For the hyperbolic target location method, target position is estimated by 
minimizing the sum of squares of difference between time delay measurements and the 
estimated function. The range equivalent time delay measurement is defined by Eq. 12 
and rewritten as 

ckT rt - (Rj. +R r )-L. (65) 

For a transmitter-target-receiver pair, baseline L is constant; therefore, the time 
delay measurement will depend on the range sum (Rj.+R r ). As in Eq. 15, range sum 

equals twice the semi-major axis length. From the definition of an ellipse, the iso-contour 
of the constant time delay measurement is elliptical with the transmitter and receiver 
position as ellipse foci (Figure 45). 
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Figure 45. Elliptical iso-contours of constant time delay measurements with 

transmitter and receiver at ellipse foci. 


In the hyperbolic target location method, an erroneous time delay measurement 
from a transmitter-target-receiver pair defines an ellipse with “thickness” in which the 
“thickness” is determined by the amount of error in the time delay measurement. The 
intersection of all ellipses approximates the uncertainty associated with the target position. 
For geometry (1), the target position estimates from the hyperbolic method result in a 
relatively circular and small error ellipse (Figure 46), while for geometry (2), the error 
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ellipse is eccentric with the semi-major axis perpendicular to the line extending from the 
target towards the receiver (Figure 47). 



Figure 46. Uncertainty area associated with position estimates from hyperbolic 

method for target in Tx-Rx cluster. 



Figure 47. Uncertainty area associated with position estimates from hyperbolic 

method for target outside Tx-Rx cluster. 
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Comparing the scatter plots from the elliptical and hyperbolic methods in Figure 
41 and Figure 42, respectively, it is evident that uncertainty associated with position 
estimates from the hyperbolic method is smaller. The error distribution of position 
estimates from elliptical and hyperbolic methods appear to be orthogonal and are more 
pronounced in geometries that exhibit DOP; that is, the semi-major axis of the error 
ellipse of estimates from the elliptical method is approximately aligned with the semi¬ 
minor axis of the error ellipse of estimates from the hyperbolic method, and vice versa. A 
possible option to improve target position estimates is to combine the elliptical and 
hyperbolic method in such a way that capitalizes on precision in their respective semi¬ 
minor axis direction. 

Every target position estimate is accompanied by a measure of confidence. The 
method of computing the error ellipse parameters outlined in Section II.E and Appendix 
A is an indication of the estimate’s accuracy. The error ellipse shows how measurement 
errors translate to the estimate’s uncertainty given certain geometry with the ellipse semi¬ 
major and semi-minor axes being directions with greatest and least deviation. The size of 
the ellipse is determined by the critical chi-square value xl, a suc ^ that the area to the left 
(l-a)xl00% is its corresponding confidence level. For a 90% confidence error ellipse, 
the critical chi-square value is 4.60, and it will enclose approximately 90% of the Monte 
Carlo position estimates. For the transmitters-target-receiver geometry in Figure 42, the 
error ellipse at various confidence levels for hyperbolic position estimates is shown in 
Figure 48. It was verified that the number of Monte Carlo estimates that fall within the 
ellipse agrees with the confidence level of the ellipse. 
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Figure 48. Hyperbolic target location estimate’s error ellipse at 50%, 70%, 90%, 

and 99% confidence levels. 
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V. CONCLUSION 


A. SUMMARY OF FINDINGS 

This thesis explored two areas pertaining to the use of multiple passive non- 
cooperative opportunistic transmitters in target detection and location. To investigate 
these PBR issues, detection coverage and target location estimation models have been 
developed using FEKO and MATLAB software. The models were used to simulate 
detection performance and location estimation accuracy for various transmitter-target- 
receiver configurations in a maritime scenario with a low-RCS target. 

The proposed PBR approach offers potential advantage in the detection of 
stealthy, low-RCS targets which are designed to minimize monostatic radar echoes. RCS 
returns of stealthy ships can be sufficiently large at certain return directions, allowing 
PBR to capitalize on this characteristic to detect low-RCS targets using opportunistic 
transmissions. Being a passive system allows the PBR receiver to remain covert, making 
it more resilient to detection and EA. The ability to leverage available transmission and to 
detect targets passively also serves to enhance situational awareness, thereby advancing 
one’s position in the OODA loop during an operation. It is also advantageous to use 
multiple transmitters at separate locations for detection and target location estimation as it 
adds spatial diversity and reduces the effect of DOP. 

1. Findings for Detection Coverage 

Target detection performance varies with the number of transmitters, their 
parameters, and receiver parameters. It also depends upon the bistatic RCS scattering 
characteristic of the target. Detection performance—received target echo SNR and area 
with good target echo SNR—improves with the number of opportunistic transmitters and 
is observed to vary significantly with a slight change in receive angle due to the rapid 
variation in RCS returns versus angle from a target. Overall detection coverage area 
increases as transmitter(s) to receiver range increase at the expense of overall sensitivity. 
Detection coverage is generally the best when the target is located within the 
transmitters-receiver cluster with the highest SNR observed when the transmitter-target- 
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receiver is in the forward-scattering configuration. Detection gaps—defined as areas with 
the low SNR at the receiver—occur when the target is outside the transmitters-receiver 
cluster or away from bearings extending from the transmitter-receiver baseline. Outside 
the cluster of transmitters-receiver, detection is better when the transmitter-target-receiver 
is in the back-scattering configuration. Other than the areas that satisfy a back-scattering 
configuration, it is observed that there are additional areas or hyperbolic-like curves— 
typical of TDOA measurements—extending from the transmitters that exhibit good 
detection coverage. Therefore, the optimal transmitters-target-receiver configuration and 
geometry for target detection coverage is such that the target is within the spatial “net” of 
transmitters-receiver. The findings with regard to detection coverage align with the 
findings and observations in [4, 40]. 

2. Findings for Target Location 

Two target location estimation methods were investigated: (1) elliptical method 
and (2) hyperbolic method. The elliptical method solves for the bistatic triangle 
parameters for each PBR pair; these parameters are then used to define transmitter-to- 
target and receiver-to-target bearings, after which the “best fit” point of intersection is 
estimated as the target position. In the hyperbolic method, time delay measurements from 
all PBR pairs are used to estimate the target’s location by finding the least-squares 
solution. The results from the target location model showed that the DF measurement 
error is the dominant source of error for the elliptical method, while the time delay 
measurement error is the primary source of error for the hyperbolic method. The 
hyperbolic method generally estimates the target’s location more accurately than the 
elliptical method as it uses fewer measurement sources. It was also observed that overall 
location accuracy is poor when measurement errors are modeled as dependent on the 
SNR of the received target echo. 

When SNR-dependent measurement errors are considered, results from the target 
location model agree closely with the observations made from the detection coverage 
results. Target location methods produce the most accurate results when the target lies 
within the transmitters-receiver cluster in a receiver-centered geometry. The transmitters- 
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receiver cluster area can be expanded by increasing the transmitters to receiver range; 
however, this reduces the overall location accuracy within the cluster. 

Monte Carlo simulations are performed for both target location methods to 
explore the relation between spatial error distribution and transmitter-target-receiver 
geometry. The results from Monte Carlo simulations in a scatter plot reveal the 
uncertainty area associated with the target location estimate for a particular geometry. 
Across all geometries, the uncertainty area associated with estimates from the hyperbolic 
method covers a smaller area compared with that from the elliptical method, suggesting 
that the hyperbolic method estimates the target location with higher accuracy. For targets 
within the cluster of transmitters-receiver, the uncertainty area would be approximately 
circular. However, when targets are outside the cluster of transmitters-receiver, the 
uncertainty area is elliptical. Graphical illustrations were presented to understand how 
measurement errors translate to uncertainty estimates for both methods. Monte Carlo 
simulation results from the target location model showed that the elliptical and hyperbolic 
target location methods have approximately orthogonal error distribution, which are more 
pronounced when the uncertainty area associated with the estimates are elliptical. That is, 
the semi-major axis of the error ellipse of estimates from the elliptical method is 
approximately aligned with the semi-minor axis of the error ellipse of estimates from the 
hyperbolic method, and vice versa. 

The method to compute the error ellipse parameters associated with the 
hyperbolic method was presented and included in the target location model. Every target 
position estimated from the hyperbolic method would be accompanied by a measure of 
confidence in the form of an error ellipse. Like Monte Carlo simulations, the error ellipse 
shows how measurement errors translate to an estimate’s uncertainty. For an error ellipse 
at (l-a)xl00% confidence, (\-a) of the 100 Monte Carlo simulated estimates will, on 
average, fall within the computed ellipse. 

B. FUTURE WORK 

Detection coverage models and target location models developed in this thesis can 
be further improved to include bistatic beam-to-beam issues and complex sources of 
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measurement and data error such as sea clutter and errors associated with sea-borne 
vessel’s motion (e.g., pitch and roll). Real data can be used in the models developed in 
this study to extend the findings regarding target detection coverage and location 
accuracy. An extension to the current work is to investigate the feasibility and 
performance of using multiple passive cooperative opportunistic transmitters in maritime 
target detection and location as proposed in [40]. A combined target location method that 
leverages on the spatial precision of the elliptical and hyperbolic target location methods 
can be explored to improve estimation accuracy. Alternatively, the use of other sources of 
measurements, like Doppler information to estimate the target’s location, can be 
investigated. 

With target detection and location in place, the next step is to explore tracking 
algorithms, classification methods, and identification methods. The MIMO tracking 
algorithm, which uses the coherent-on-receiver method [41] or the track-before-detect 
method that enables tracking and detection of targets in low SNR environments [42], are 
feasible options. As stated in [4], an interesting extension would be to consider the use of 
bistatic ISAR signatures in the forward-scattering and back-scattering configurations to 
aid in target identification and classification [25]. 
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APPENDIX A. DERIVATION OF ERROR ELLIPSE PARAMETERS 
FROM BIVARIATE NORMAL DISTRIBUTION 


This appendix presents the derivation of parameters that define the error ellipse 
[38]. To understand the how measurement errors translate to estimate position uncertainty 
in the form of an error ellipse, consider a univariate normal distribution with mean // and 

variance a 1 for the random variable X : 


/(*) = 


1 

sfl ~Ka 2 


e 


(x-Mf 

2a 2 


( 66 ) 


for —oo < x < go , —oo < ju< go and cr 2 > 0. 

Similarly, a multivariate normal distribution with mean p and covariance X for 
random variables in X is expressed as 


/(X) = 


1 

J(2x) p \I\ e 


(67) 


for —OO < X < oo, —oo < p < x and for positive definite X . Eq. 67 is derived from Eq. 66 


by rewriting the term in linear algebra formulation: 

(6g) 

(x-^) r I-'(x-^). (69) 

Eq. 69 is the squared Mahalanobis distance where X and p are (pxl) vectors, 


Z is a (pxp) positive definite covariance matrix, and p denotes the number of 
variables in vector X. 

To derive the semi-major axis and semi-minor axis lengths of an error ellipse in 
the XY-plane, reduce Eq. 67 into a bivariate normal distribution where 


(70) 


x = 

X 



,1 = 

<T 2 POO 

x r xy x y 

2 




_/A_ 


poo o 

• yx y x y 


By substituting Eq. 70 into the Mahalanobis distance in Eq. 69, 
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(x-n) s- , (x-n)= — 


p 


xy 


(,x-p x f 2 p^x-ju x \x-ju y ) ( x-ju y ) 


+ - 


(7 <J 

x y 


(71) 


the probability density function of a pair of jointly Gaussian random variables is then expressed 
as 


fxY y ) 


(_ I 

(x-jU x f 2pxy ( X-jU x )(x-/i y ) 1 (x-jU y f 

[ 2(1-4) 

L ^ J 


for —oo < x < oo , —oo < y < oo. 


(72) 


Figure 49 shows the joint Gaussian probability density function (pdf) for various 
correlation and standard deviation values. The shape and orientation of the joint Gaussian 
pdf surface depends on the values of a x , <j y , and p xy , while the size of the surface 

contours depends on the pdf value f XY (x, y) in consideration. 


An error ellipse is defined as a locus of points with constant f XY (x, y) values, as 
seen in the surface contour plots in Figure 49. Graphically, the error ellipse is the points 
of intersection between the pdf surface and a horizontal plane. Given the values of <j x , 

<j , and p xy , fxY ( x > >’) is constant for values of x and y when the exponential argument 
in Eq. 72 is constant. Mathematically, this is expressed as 

(X-p) T Z' 1 (X-p) = c. (73) 

Let Xj = X- p such that Eq. Error! Reference source not found, reduces to 

X[z _1 Xj = c. 

Factorize the covariance matrix Z by eigen-decomposition: 


(74) 


Z = PDP , 


(75) 


where P is a matrix of eigenvectors and D is a diagonal matrix whose elements are the 
corresponding eigenvalues (/l l ,A 2 ). 
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Figure 49. Joint Gaussian pdf surface and contours for various <j x , o y , and p xy 

values. After [38]. 
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The eigenvectors of the covariance matrix represent the direction in which the 
data varies the most and are the vectors along the error ellipse’s semi-major and semi¬ 
minor axis. The corresponding eigenvalues indicate the spread of the data in the direction 
of the eigenvectors (Figure 50). For an axis-aligned error ellipse such that covariance is 
zero, the eigenvalues equal the variances of the covariance matrix, and the eigenvectors 
are in the direction of the x-axis and y-axis. However, for correlated data, the 
eigenvectors represent the direction with the largest spread of data and the eigenvalue 
determines the extend of spread. 



Figure 50. Eigenvectors (^A l ,A 7 ) of a covariance matrix on an error ellipse. 

If X is a symmetric matrix, P _1 = P r , Eq. 75 is then expressed as 

Z = PDP r . (76) 


When X can be eigen-decomposed and when none of its eigenvalues are zero, then X is 
non-singular and its inverse is given by 

S =PDP r (77) 

and 

D 1 = P T X 1 P. (78 ) 
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As illustrated in Figure 50, define the r-axis and .s-axis as the principle axis of the 
error ellipse. The principle axis is defined as rotating the original x-axis and y-axis 
through an angle 6 and is expressed mathematically as 


= P 


(79) 


To determine the semi-major and semi-minor axis lengths of the error ellipse, 
substitute Eq. 78 and Eq. 79 into Eq. 74 

[r 


f ] =c 
[r s]D~ l [r sf =c 


[r s] 


1 / 


% 

0 


0 


a 


[r sf=c 


■ + - 


= c. 


(80) 


\ /l 2 

The general form of an ellipse centered at the origin is 

2 2 

^ + ^ = 1 

a 2 b 2 ( 81 ) 

where a and b are the semi-major and semi-minor axis lengths, respectively. Comparing 
Eq. 80 and Eq. 81, 

a = ^c\\ 

I( 82 ) 

To understand the significance of the constant c , consider a (px 1) multivariate 
normal vector X ~ N(\i, X) . Geometrically, degrees of freedom is interpreted as the 
dimension of the vector subspace. Hence, the squared Mahalanobis distance between X 
and n is a chi-square distribution with p degrees of freedom (Figure 51): 

\ T V-l , 


(X-|») Z-^x-ii)-^. 


(83) 
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Chi-Square PDF with p Degrees of Freedom 



Figure 51. Chi-square pdf for p degrees of freedom. 


For a hyper-ellipse, defined as the squared Mahalanobis distance equal to a 
critical value of a chi-square distribution with p degrees of freedom and evaluated at a , 
the probability of the random vector X falling inside the ellipse is 1 — a . 


Prj(X- t i) r I- | (X- (l )<^[ = l- 


■a 


(84) 


where a refers to the area to the right of the chi-square critical value (Figure 52). 

Therefore, the hyper-ellipse is defined as a locus of points with constant 
fxY (x, y) values such that 

(x-rt r z- 1 (x-rt = ^. 

By comparing Eq. 73 and Eq. 85, Eq. 82 can be written as 

a = \[x2^a\ 


(85) 




'( 86 ) 
For a two-dimensional hyper-ellipse, the chi-squared distribution with two 
degrees of freedom, p = 2 will be used. 
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Chi-Square PDF with 2 Degrees of Freedom 



Figure 52. Chi-square pdf where the area to the right of the critical value is a . 

For an error ellipse with confidence level (1 -a)x\ 00%, the relevant x\, a va ' ue 

in Eq. 86 is obtained from the chi-square distribution table. A sample of a chi-square 
distribution table is given in Table 8. 


Table 8. Chi-square distribution table 


Degrees of 
Freedom 

Probability 

0.95 

0.90 

0.80 

0.70 

0.50 

0.30 

0.20 

0.10 

0.05 

0.01 

0.001 

1 

0.004 

0.02 

0.06 

0.15 

0.46 

1.07 

1.64 

2.71 

3.84 

6.64 

10.83 

2 

0.10 

0.21 

0.45 

0.71 

1.39 

2.41 

3.22 

4.60 

5.99 

9.21 

13.82 

3 

0.35 

0.58 

1.01 

1.42 

2.37 

3.66 

4.64 

6.25 

7.82 

11.34 

16.27 

4 

0.71 

1.06 

1.65 

2.20 

3.36 

4.88 

5.99 

7.78 

9.49 

13.28 

18.47 

5 

1.14 

1.61 

2.34 

3.00 

4.35 

6.06 

7.29 

9.24 

11.07 

15.09 

20.52 


Referring to Figure 53, the error ellipse’s principle axis [r s] is defined as the 

counterclockwise rotation of the original xy-axis [x y] through an angle 6. The vector 

representation of the original xy-coordinates as a function of the ^-coordinates is 

q = Ay (87) 
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where c, = [r 



y] r and A 


cos # sin # 
-sin# cos# 



Figure 53. Rotation of error ellipse principle axis. From [38], 


From Eq. 87, 

=A I y A r , 

y O' 

? L° <’ 

cr 2 poo 

Z x r'xy x y 

y ~ 2 

P O O O 

_r xy x y y 

Substituting Eq. 89 and Eq. 90 into Eq. 88 yields 


o 

_1 


cos# 

sin# 


poo 

r'xy x y 

cos# 

-sin# 

1- 

o 

a 

1_ 


-sin# 

cos# 

Pxy^x^y 

J 

sin# 

cos# 


( 88 ) 

(89) 

(90) 
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which is further reduced to 


OOO 00 

cr. = a x cos 0+2.p a x a sin 6 cos6+<j y sin 6 , 

0 0 0 0 0 

cr = (j v sin 0-2cr cr sin 6 cos 6- cr cos 0, 

j a " xy jc y y 

0 = (cr 2 - cr j sin 0cos 0 + p xy <r x <j y (cos 2 0-sin 2 0). 


Simplifying Eq. 92 gives the orientation of the error ellipse: 


0 = -tan 1 
2 


2 P, 


cr cr 

x e y e x e y e 


V 


(91) 

(92) 

(93) 


(94) 


85 



THIS PAGE INTENTIONALLY LEFT BLANK 


86 



APPENDIX B. SPECIFICATION SHEET FOR MANTADIGITAL 
RADAR BY KELVIN HUGHES 


NAVIGATION 
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MANTADIGITAL™ RANGE 



Figure 54. Technical Specifications for MantaDigital Radar by Kelvin Hughes. 

From [43], 
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APPENDIX C. SPECIFICATION SHEET FOR TELEDYNE 
DEFENCE QR026 EW RECEIVER 
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Figure 55. Technical Specifications for Teledyne Defence QR026 EW Receiver. 

From [44], 
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APPENDIX D. SPECIFICATION SHEET FOR POYNTING 
DEFENCE DF-A0062 DF RECEIVER 




Poynting 

Antennas 


VERSION: 2.0 


Direction Finding and 
Monitoring Antenna 

20 - 6000 MHz _ 

Product Code: DF-A0062 



PRODUCT DESCRIPTION 

Th*s direction finding antenna cover* a frequency range 
of 20 Mdz lo 6 GHz. A separate. but fully ntegraled 
active mentoring anaenrta gves high sansitrvty onrvn- 
dineoccnai coverage on the same a us at the DF 
anfienna without nterterng .ith the DF array 

The tuli-size dement* on all band* give excellent DF 
sensrtr.ity jttmare angular resolution for strong sigrars 
is wel under 1* for most of the frequency range Diode 
dements protide good cross-pci ansa bon rejection and 
fair performance tor srgras am*ng *cm up to I5f* 
above or betov. me horizon 

The ntegated montonnp antenna is In two bands, 
mounted above me Bard C S D OF antenna nrxxJUe 
boused in the raceme Each band of the montonng 
array ts amplified at the top of me cabe and passr.e > 
oombmed to give continuous coverage o%e-r the 
frequency range 20 fifrfz to 6 GHz. wflh a single cutout 
connector 


SPECIFIC ATKJNS 


Product codes: 


DFAJC062 

DF v.ith monhoma up to 6 GHz 



Electrical OF 


DF frequency range 

Band A 20 - 300 MHz; 

Band B 300-1000 MHz 

EandC iCOa - >300 MHz. 

Eana D. 3000 - 6300 MHz 

Normal nput impedance 

r 50Q 

Antenna type 

5-etement DF interferometer 
i From 5 to 2-channel reoerver 
ccmpattbiei 

Polarisation 

Vertical 

DjtpU : 

RG 400 cables 1 

colors 

20 x TNC male 

Mowtorog connector 

1 »N ma* 



B>Ctficai - monitoring 


Frequency range 

23 - 6000 MHz 

Norma* H-piane beam.v»dth 

360' 

Normal £ -pane beamwidm 


Typical '/SWR* 

2t 

Potansabon 

Vertical 

Connector type 

N-type male 


' W ' 1 

Ncmnai mpecance 

r >p‘* .v.Tajv , ca v.a> 

13-24 V DC 

input current* 

< tiO mA 

Po- er ccnsu-npOcn ;ncmr 

rt 225 W 

CPIdB itvo t* 

>It dBm 

OP 2 ihpjr__ 

> 31 dBm 

IWESp 

» 21 dBm 

Sensitruty 

<typi 

\SN* 

ooe sw« 

1 Kzi* 

20- 100 MHz 

-20 dBpY.yn 

1^1000 

MHz 

•30 dBpV.tn 

1 - 3 GHz 

-35 dBuVihn 

3 - 6 GHZ 

-40 dBpv.hi 



Mechanical: 

Crosvsectiqi 
Maxmum *+ 

nai mM area 

• 

r>(j $i>eed 

160 kmh iwehout *ce load) 

Assembled height 

[>60 m 

Assembled dameter impi 

2 70 m 

Snapping dmensons 
iVeighr ort artenna 
ncsudng shipping container 

<V) kg 

135 kg 


Tha DF amenra is designed to be usable ,tth either a 
5- or 2-channel phese-sererthe receiver and corteiaflve 
aigcrdhm Char ec!ensaeon of the antenna can be 
performed an request 

Related products DF-A0064 CF-A0357-03. DF-AO038 
and DF-AOOS6 


ELECTRICAL FEATURES: 

• Ftihsize DF 

• 'lYd* bard OF 

• S-eiemerr iraerforcmatar 

• Cptmaed for 2-channel receivers 

• High sensitive, omm artenna integrated 

• Tactical rapid deploy sto.v design .-tm integrated oust caps 


PAGE 1 of 4 



MECHANICAL FEATURES: 

• Designed tor tov,«r maurftng 

defencegpoynbng oo z& 

www.poyntingdefence.com 

► — u » fiMfi i 


DEFENCE and 
SPECIALISED 


Figure 56. Technical Specifications for Poynting Defence DF A0062 DF 

Receiver. From [45], 
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APPENDIX E. DETECTION COVERAGE PLOTS FOR TARGET 

PLANE IN S-BAND 


Referring to detection coverage for the plate at 0° and 180° orientations in Figure 
57 and Figure 58, respectively, received SNR is the greatest when the transmission and 
reception of the target pane echo occurs in the vicinity of the target plane faces (at 90° 
and 270° RCS mainlobes) and minimum when transmission and reception occurs at the 
edges of the plate (0° and 180°). Similar observations were made for detection coverage 
for target plane at 90° and -90° orientation (Figure 59 and Figure 60). 

Comparing the detection coverage for a plate at a 45° orientation with different 
transmitter positions in Figure 61 and Figure 62, it is evident that different transmitter- 
receiver geometry can result in similar detection coverage. The target position with the 
greatest SNR is observed when either transmission or reception LOS falls in the vicinity 
of the plate faces. Similar observations were made for a plate at a -45° orientation (Figure 
63 and Figure 64). 



Figure 57. Detection coverage (left) for Tx at [-5000, 0] and target plate at 0° 
orientation with corresponding S-band bistatic RCS (right). 
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figure 58. 


Detection coverage (left) for Tx at [-5000, 0] and target plate at 0° 
orientation with corresponding S-band bistatic RCS (right). 


Detection Coverage SNR (dB) 



Figure 59. Detection coverage (left) for Tx at [-5000, 0] and target plate at 90° 
orientation with corresponding S-band bistatic RCS (right). 
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Figure 60. Detection coverage (left) for Tx at [-5000, 0] and target plate at -90° 
orientation with corresponding S-band bistatic RCS (right). 
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Figure 61. Detection coverage (left) for Tx at [-5000, 0] and target plate at 45° 
orientation with corresponding S-band bistatic RCS (right). 


95 





Distance (m) 


Detection Coverage SNR (dB) 


4 


x 10 



Figure 62. Detection coverage (left) for Tx at [0, 5000] and target plate at 45° 
orientation with corresponding S-band bistatic RCS (right). 




Figure 63. Detection coverage (left) for Tx at [-5000,0] and target plate at -45° 
orientation with corresponding S-band bistatic RCS (right). 
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Figure 64. Detection coverage (left) for Tx at [0, 5000] and target plane at -45° 
orientation with corresponding S-band bistatic RCS (right). 
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APPENDIX F. DETECTION COVERAGE RESULTS 


F.l VARY NUMBER OF TRANSMITTERS 



*2 '15 -1 -05 0 05 1 15 2 

Ostanoe <m) , 


Figure 65. Detection coverage at S-Band for 1 Tx at 5 km range and target at 0° 

orientation. 



CHtanoe <m) , jq* 


Figure 66. Detection coverage at S-Band for 1 Tx at 5 km range and target at 

-90° orientation. 
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L>-starve* <m) , 1 q* 

Figure 67. Detection coverage at X-Band for 1 Tx at 5 km range and target at 0° 

orientation. 



CHtanoe <m) x y q 

Figure 68. Detection coverage at X-Band for 1 Tx at 5 km range and target at 

-90° orientation. 
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Figure 69. Detection coverage at S-Band for 2 Tx at 5 km range and target at 0° 

orientation. 



CVstanoe <m> , j q 


Figure 70. Detection coverage at S-Band for 2 Tx at 5 km range and target at 

-90° orientation. 


101 



-2 '15 -1 -05 0 05 1 15 2 

CHtanoe <m) , iq* 


Figure 71. Detection coverage at X-Band for 2 Tx at 5 km range and target at 0° 

orientation. 



Figure 72. Detection coverage at X-Band for 2 Tx at 5 km range and target at 

-90° orientation. 
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Ostanoe <m> , 

Figure 73. Detection coverage at S-Band for 3 Tx at 5 km range and target at 0° 

orientation. 



Orslanoe <m) , 1 q < 


Figure 74. Detection coverage at S-Band for 3 Tx at 5 km range and target at 

-90° orientation. 
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Figure 75. Detection coverage at X-Band for 3 Tx at 5 km range and target at 0° 

orientation. 



CHtanoe <m) , f q* 


Figure 76. Detection coverage at X-Band for 3 Tx at 5 km range and target at 

-90° orientation. 
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1)-starve* <m) , 


Figure 77. Detection coverage at S-Band for 4 Tx at 5 km range and target at 0° 

orientation. 



CHtanoe <m) , f q* 


Figure 78. Detection coverage at S-Band for 4 Tx at 5 km range and target at 

-90° orientation. 
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Detedton Coverage SNR <dB> 
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Figure 79. Detection coverage at X-Band for 4 Tx at 5 km range and target at 0° 

orientation. 



CHtanoe <m) K y q 

Figure 80. Detection coverage at X-Band for 4 Tx at 5 km range and target at 

-90° orientation. 
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Figure 81. Detection coverage at S-Band for 8 Tx at 5 km range and target at 0° 

orientation. 



CHtanoe <m) , 1 q* 


Figure 82. Detection coverage at S-Band for 8 Tx at 5 km range and target at 

-90° orientation. 
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Detection Coverage SNR <dB> 
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Figure 83. Detection coverage at X-Band for 8 Tx at 5 km range and target at 0° 

orientation. 
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Figure 84. Detection coverage at X-Band for 8 Tx at 5 km range and target at 

-90° orientation. 
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F.2 VARY TRANSMITTER RANGE 



Balance im> x ^ q * 


Figure 85. Detection coverage at S-Band for 1 Tx at 5 km range and target at 0° 

orientation. 


Detection Coverage SNR (dBl 
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Figure 86. Detection coverage at S-Band for 1 Tx at 15 km range and target at 0° 

orientation. 
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Detection Coverage SNR (dB| 



Figure 87. Detection coverage at S-Band for 2 Tx at 5 km range and target at 0° 

orientation. 


Detection Coverage SNR (dB| 
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Figure 88. Detection coverage at S-Band for 2 Tx at 15 km range and target at 0° 

orientation. 
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Detection Coverage SNR (dB| 



Figure 89. Detection coverage at S-Band for 4 Tx at 5 km range and target at 0° 

orientation. 
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Figure 90. Detection coverage at S-Band for 4 Tx at 15 km range and target at 

0° orientation. 
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Figure 91. Detection coverage at S-Band for 8 Tx at 5 km range and target at 0° 

orientation. 



Figure 92. Detection coverage at S-Band for 8 Tx at 15 km range and target at 0° 

orientation. 
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F.3 VARY TRANSMITTER-TARGET-RECEIVER GEOMETRY 



Figure 93. Detection coverage at S-Band for 4 Tx in receiver-centered geometry 

and target at 0° orientation. 


Detection Coverage SNR (dB| 



Figure 94. Detection coverage at S-Band for 4 Tx in receiver-centered geometry 

and target at -90° orientation. 
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Detection Coverage SNR (dB) 



Figure 95. Detection coverage at S-Band for 4 Tx in receiver-centered geometry 

and target at 45° orientation. 



Figure 96. Detection coverage at S-Band for 4 Tx in receiver-centered geometry 

and target at 30° orientation. 
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Figure 97. Detection coverage at S-Band for 4 Tx in transmitter-clustered 

geometry and target at 0° orientation. 



Figure 98. Detection coverage at S-Band for 4 Tx in transmitter-clustered 
geometry and target at -90° orientation. 
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Figure 99. Detection coverage at S-Band for 4 Tx in transmitter-clustered 

geometry and target at 45° orientation. 
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Figure 100. Detection coverage at S-Band for 4 Tx in transmitter-clustered 

geometry and target at 30° orientation. 
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F.4 RANDOMLY DISTRIBUTED TRANSMITTERS 



[>*tartoe |ftl> x 5Q 4 

Figure 101. Detection coverage at S-Band for 8 Tx randomly positioned around 
Rx at 5-20 km range and target at 0° orientation. 


Detection Coverage SNR (dB| 



Figure 102. Detection coverage at S-Band for 8 Tx randomly positioned around 
Rx at 5-30 km range and target at 0° orientation. 
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3, iq 4 Detection Coverage SNR «iB> 



Figure 103. Detection coverage at S-Band for 20 Tx randomly positioned around 
Rx at 5-20 km range and target at 0° orientation. 


t io 4 Detection Coverage SNR «JB> 



Figure 104. Detection coverage at S-Band for 20 Tx randomly positioned around 
Rx at 5-30 km range and target at 0° orientation. 
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Detection Coverage SNR (dBi 



Figure 105. Detection coverage at S-Band for 8 Tx randomly positioned around 
Rx at 5-20 km range and target at -90° orientation. 



Distance im> 


Figure 106. Detection coverage at S-Band for 8 Tx randomly positioned around 
Rx at 5-30 km range and target at -90° orientation. 
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Figure 107. Detection coverage at S-Band for 20 Tx randomly positioned around 
Rx at 5-20 km range and target at -90° orientation. 


3, jo* Detection Coverage SNR «JB> 



Figure 108. Detection coverage at S-Band for 20 Tx randomly positioned around 
Rx at 5-30 km range and target at -90° orientation. 
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Figure 109. Detection coverage at S-Band for 8 Tx randomly positioned around 
Rx at 5-20 km range and target at 45° orientation. 



Figure 110. Detection coverage at S-Band for 8 Tx randomly positioned around 
Rx at 5-30 km range and target at 45° orientation. 


121 







Detection Coverage SNR <dB> 



Figure 111. Detection coverage at S-Band for 20 Tx randomly positioned around 
Rx at 5-20 km range and target at 45° orientation. 



Figure 112. Detection coverage at S-Band for 20 Tx randomly positioned around 
Rx at 5-30 km range and target at 45° orientation. 
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APPENDIX G. TARGET PATH DETECTION RESULTS 



SNR(dB) 



10 
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Figure 113. Detection coverage along target’s path (2 km resolution) at S-Band 

for 4 Tx. RCS computed by calling FEKO. 


SNR at Rx along Frigate Path 



Figure 114. SNR at receiver along target’s path (2 km resolution) at S-Band for 4 

Tx. RCS computed by calling FEKO. 
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Figure 115. Detection coverage along target’s path (100 m resolution) at S-Band 
for 4 Tx. RCS extracted from pre-computed RCS table. 


SNR at Rx along Frigate Path 



Figure 116. SNR at receiver along target’s path (100 m resolution) at S-Band for 
4 Tx. RCS extracted from pre-computed RCS table. 
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Figure 117. Detection coverage along target’s path (100 m resolution) at S-Band 
for 8 Tx randomly position. RCS extracted from pre-computed RCS 

table. 


SNR at Rx along Frigate Path 



Figure 118. SNR at receiver along target’s path (100 m resolution) at S-Band for 
8 Tx randomly position. RCS extracted from pre-computed RCS table. 
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APPENDIX H. TARGET LOCATION RESULTS 
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Figure 119. Target location errors using elliptical method (S-Band, 4 Tx at 5 km 
range) and SNR-independent measurement errors. 
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Figure 120. Target location errors using hyperbolic method (S-Band, 4 Tx at 5 
km range) and SNR-independent measurement errors. 
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Figure 121. Target location errors using elliptical method (S-Band, 8 Tx at 5 km 
range) and SNR-independent measurement errors. 
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Figure 122. Target location errors using hyperbolic method (S-Band, 8 Tx at 5 
km range) and SNR-independent measurement errors. 
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Figure 123. Target location errors using elliptical method (S-Band, 4 Tx at 15 
km range) and SNR-independent measurement errors. 
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Figure 124. Target location errors using hyperbolic method (S-Band, 4 Tx at 15 
km range) and SNR-independent measurement errors. 
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Figure 125. Target location errors using elliptical method (S-Band, 8 Tx at 15 
km range) and SNR-independent measurement errors. 
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Figure 126. Target location errors using hyperbolic method (S-Band, 8 Tx at 15 
km range) and SNR-independent measurement errors. 
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Tarpet Pos*»on Estimate Error (m) 
using Tn geometry and Least-Squares Method 
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Figure 127. Target location errors using elliptical method (S-Band, 4 Tx 
clustered) and SNR-independent measurement errors. 
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Figure 128. Target location errors using hyperbolic method (S-Band, 4 Tx 
clustered) and SNR-independent measurement errors. 
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Figure 129. Target location errors using elliptical method (S-Band, 8 Tx 
clustered) and SNR-independent measurement errors. 
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Figure 130. Target location errors using hyperbolic method (S-Band, 8 Tx 
clustered) and SNR-independent measurement errors. 
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Target Position Estimate Error (m) 



Figure 131. Target location errors using elliptical method (S-Band, 4 Tx at 5 km 
range) and SNR-dependent measurement errors. 
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Figure 132. Target location errors using hyperbolic method (S-Band, 4 Tx at 5 
km range) and SNR-dependent measurement errors. 
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Figure 133. Target location errors using elliptical method (S-Band, 8 Tx at 5 km 
range) and SNR-dependent measurement errors. 
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Figure 134. Target location errors using hyperbolic method (S-Band, 8 Tx at 5 
km range) and SNR-dependent measurement errors. 
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Figure 135. Target location errors using elliptical method (S-Band, 4 Tx at 15 
km range) and SNR-dependent measurement errors. 
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Figure 136. Target location errors using hyperbolic method (S-Band, 4 Tx at 15 
km range) and SNR-dependent measurement errors. 
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Figure 137. Target location errors using elliptical method (S-Band, 8 Tx at 15 
km range) and SNR-dependent measurement errors. 
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Figure 138. Target location errors using hyperbolic method (S-Band, 8 Tx at 15 
km range) and SNR-dependent measurement errors. 
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Figure 139. Target location errors using elliptical method (S-Band, 8 Tx at 
random positions) and SNR-dependent measurement errors. 
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Figure 140. Target location errors using hyperbolic method (S-Band, 8 Tx at 
random positions) and SNR-dependent measurement errors. 
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target Position Estimate Error (m) 
using Tngormetry and Least-Squares Method 
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Figure 141. Target location errors using elliptical method (S-Band, 4 Tx 
clustered) and SNR-dependent measurement errors. 
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Figure 142. Target location errors using hyperbolic method (S-Band, 4 Tx 
clustered) and SNR-dependent measurement errors. 
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Figure 143. Target location errors using elliptical method (S-Band, 8 Tx 
clustered) and SNR-dependent measurement errors. 
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Figure 144. Target location errors using hyperbolic method (S-Band, 8 Tx 
clustered) and SNR-dependent measurement errors. 
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